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ABSTRACT 

We show that it is possible to improve the chiral behaviour and the 
approach to the continuum limit of correlation functions in lattice 
QCD with Wilson fermions by taking arithmetic averages of corre- 
lators computed in theories regularized with Wilson terms of opposite 
sign. Improved hadronic masses and matrix elements can be obtained 
by similarly averaging the corresponding physical quantities separately 
computed within the two regularizations. To deal with the problems 
related to the spectrum of the Wilson-Dirac operator, which are par- 
ticularly worrisome when Wilson and mass terms are such as to give 
contributions of opposite sign to the real part of the eigenvalues, we 
propose to use twisted-mass lattice QCD for the actual computation 
of the quantities taking part to the averages. The choice ±7r/2 for the 
twisting angle is particularly interesting, as 0(a) improved estimates of 
physical quantities can be obtained even without averaging data from 
lattice formulations with opposite Wilson terms. In all cases little or 
no extra computing power is necessary, compared to simulations with 
standard Wilson fermions or twisted-mass lattice QCD. 
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1 Introduction 



The two major limitations of the standard Wilson formulation of lattice 
QCD (LQCD) [Ij s-re related to the non-negligible magnitude of the system- 
atic errors associated to discretization effects and the complicated pattern 
of operator mixing induced by the explicit breaking of chiral symmetry due 
to the presence of the Wilson term in the lattice action [21 El 13 • Conse- 
quently a lot of effort has gone in the direction of reducing errors related 
to the finiteness of the lattice spacing, a, with the full implementation of 
the Symanzik program |S] up to 0(a) corrections [01 IZllHlini Ej and in de- 
veloping formulations in which chiral symmetry is exact at the level of the 
regularized action [TTIIT^OITH . 

Elimination of 0(a) effects can be in principle fully achieved. The stan- 
dard approach developed in refs. [HIIZIIHI El ^] requires, besides a suitable 
improvement of operators, the introduction in the standard Wilson LQCD 
action of the so-called Sheikholeslami-Wohlert (SW) term, 

•S'sw = a^csw(9o) Yl > (1-1) 

X 

with an appropriate value of the coefficient function csw(fl'o) depending on 
the detailed form of the pure gauge action and Wilson term ^. All the above 
modifications, particularly in the unquenched case, make 0(a) improved 
computations significantly more CPU-time demanding than standard Wil- 
son simulations. This is even more so if the evaluation of off-shell or gauge 
non- invariant correlators is attempted jlUj . 

Fully chiral invariant lattice actions have been also constructed \12\ 
I14| . where the Ginsparg- Wilson relation is exactly implemented and 
a first exploration of their numerical potentiality has been carried out |16| 
1171 118j . Overlap, domain- wall and fixed-point fermions have been put at 
work with rather encouraging results, though the extra computational effort 
required with respect to the standard Wilson case is in most cases very 
large. Of course this comparison is not totally fair because the quality of 
data produced with fermions obeying the Ginsparg-Wilson relation looks 
superior to that usually obtainable with standard Wilson fermions. 

If one does not want to make use of chirally invariant fermions, because 
of their relatively heavy computational cost, a nice alternative is to use 
twisted-mass lattice QCD (tm-LQCD) |19| l2()j. which represents a cheap but 

^In eq. 11. Ill is any lattice discretization of the gauge field strength |15| . go is the 
bare gauge coupling and we define cr^^ = i[7^,7^]/2, with {7,^,71/} = 25^^ and jI = 7^. 
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effective solution to the problem of exceptional configurations occurring in 
quenched (or partially quenched) simulations. Recently it has been shown, in 
fact, that tm-LQCD can be judiciously used to alleviate both the numerical 
difficulties associated to 0(a) discretization errors [HI and the mixing 
problem of Wilson fermions in the calculation of and the CP-conserving 
A5 = 1 weak matrix elements 

The present one is the first of two companion papers in which we pro- 
pose a simple strategy to get from simulations employing Wilson fermions 
(or rather, as we shall see, twisted Wilson fermions) lattice data that are 
free of 0(a) discretization effects and have a somewhat smoother and more 
chiral behaviour near the continuum limit than data obtained in unimproved 
standard Wilson or tm-LQCD simulations. 

In this paper we show that, if quarks are arranged in SU{2)f flavour dou- 
blets, 0(a) discretization effects are absent from the average of correlators 
(Wilson average, W A) computed with lattice actions having Wilson terms 
of opposite sign and a common value of the subtracted (unrenormalized) 
lattice quark mass 



where Mq is the bare mass. For short in the following we will occasionally 
call mq the "excess" quark mass. 

Absence of 0(a) discretization errors in 1^^'s is proved by referring 
to the Symanzik expansion of (connected, on-shell) lattice correlators in 
terms of continuum Green functions and exploiting the relations derived by 
matching the "T^s-parity" of lattice correlators under the transformation 



to the 7^5-parities of the related continuum Green functions. 7^5 is an el- 
ement of the chiral SU{2)l ® SU{2)b, group, for it can be expressed as 
the product of the following three transformations Uy^Tr) = exp(i7rri/2), 
tiy (vr) = exp(z7rT2/2), u\ {tt) = exp(i7r75r3/2), with ti,T2,ts the Pauli 
matrices. Obviously [T^s]^ = 1. 

The improvement of VF^'s is the result of the use of simple dimen- 
sional and symmetry considerations combined with certain algebraic iden- 
tities holding between pairs of correlators computed with opposite signs of 
Mq and r. These identities come from a generalized spurion analysis, which 
in turn follows from the fact that the lattice fermionic action is left invariant 
if the field transformation (|1.3j) is accompanied by a simultaneous change of 



mg = Mq- M, 



(1.2) 




(1.3) 
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sign of Mq and the coefficient, r, sitting in front of the Wilson term. Prom 
these considerations one can prove that the Symanzik coefficients necessary 
to match lattice correlators to the continuum ones have such properties un- 
der r — > — r that all 0(a) corrections to the continuum result cancel in 
WA's. 

An important point to stress is that, in order for the formal counting of 
powers of a yielded by the relevant Symanzik expansion to be meaningful, 
it is necessary that one is dealing with expectation values of multiplicatively 
renormalizable (m.r.) lattice operators. 

The obvious problem with the approach we have described is the exis- 
tence of real eigenvalues of the massless Wilson-Dirac (WD) operator with 
magnitude much smaller than the inverse lattice cutoff. For realistic values 
of the lattice spacing, the existence of such spurious eigenvalues (in partic- 
ular those that are negative at niq = 0) has been shown to be at the origin 
of the occurrence of zero modes of the massive (m^ > 0) WD operator ^ on 
the so-called "exceptional" configurations UHlllZl- The problem becomes 
more severe when the relative sign of the coefficients of Wilson and quark 
mass terms is such that they give opposite contributions to the real part 
of the spectrum of the WD operator. The remedy to this situation, which 
we will for short refer to as the problem of "spurious fermionic modes", is 
to use tm-LQCD [^|^. Fully 0(a) improved lattice data for energy lev- 
els (hence hadronic masses), matrix elements and renormalization constants 
can be obtained, without the need of computing anyone of the usual lattice 
improvement coefficients, if one takes averages of correlators evaluated with 
tm-LQCD actions having twisted Wilson terms of opposite sign. 

The procedure we propose is very flexible as it is based on the freedom of 
regularizing different flavours with Wilson terms of different chiral phases. 
Its extension to mass non-degenerate quarks is discussed in ref. [2S]. In 
the second of this series of paper 2^ we show that, without loosing 0(a) 
improvement, this freedom can be appropriately exploited to have a real and 
positive fermionic determinant, while at the same time canceling all flnite 
and power divergent contributions due to mixing with operators of "wrong 
chirality" (i.e. those due to the breaking of chiral symmetry induced by the 
presence of the Wilson term in the lattice action) in the calculation of the 
CP-conserving matrix elements of the effective weak Hamiltonian. These 
results confirm and largely extend arguments given in refs. |191l24j . 

The present paper is logically divided in two parts. In the first part 
(sections 121 and 12)) we derive our results momentarily ignoring problems with 

^For a first hint to the existence of such a problem see ref. |25|. 
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the spectrum of the WD operator. This key issue will be taken up in the 
second part (sections |1] to IHI) where, after an introduction to tm-LQCD, 
we prove that with appropriate choices of twisting none of the nice results 
derived for Wilson fermions is lost. On the contrary, if the special choice 
ib7r/2 for the twisting angle is made, simple symmetry considerations imply 
that many interesting physical quantities can be extracted from lattice data 
that are automatically 0(a) improved with no need of taking any WA ^. 

The detailed plan of the paper is as follows. In section |21 after briefly re- 
viewing the standard Wilson formulation of LQCD, we show that the critical 
mass, Mcr, changes sign if the sign in front of the Wilson term is reversed. 
This is done by providing an argument both in perturbation theory (PT) 
and in the presence of spontaneous chiral symmetry breaking (S^SB). In 
section 131 we outline the basic ideas of the present approach by showing that 
averages of lattice correlators computed in theories with Wilson terms of 
opposite sign and a common value of niq are 0(a) improved. We also prove 
that similar averages of hadronic masses, matrix elements and renormaliza- 
tion constants analogously lead to 0(a) improved quantities. We present 
in section 21 a general discussion of tm-LQCD with special attention to the 
non-singlet Ward-Takahashi identities (WTI's) of the theory. In section |21 
the idea of improvement through M^A's is extended to tm-LQCD and the 
partial breaking of parity invariance that occurs in this regularization is care- 
fully analyzed. We show that it is possible to attribute to the eigenstates of 
the transfer matrix a parity label that in the continuum limit turns out to 
coincide with the standard parity quantum number. In sect. Ill we discuss 
the remarkable simplifications that occur if the choice ib7r/2 for the twisting 
angle is made and we prove that many interesting physical quantities are 
automatically 0(a) improved. Conclusions and an outlook of future lines 
of investigations can be found in section [7J We leave in a number of ap- 
pendices certain more technical arguments and a detailed discussion of the 
improvement of the non-singlet axial WTI's for standard Wilson fermions. 

2 Wilson fermions and 75-chirality 

LQCD for one flavour doublet (A^f = 2) of quarks, as formulated by Wil- 
son is expressed in terms of the (Euclidean) partition function 

^LQCD = / V^,{^, ^, U) e-[^^^M(f/)+^P^-(^-V;,f/)] , (2.1) 

^We are indebted to M. Delia Morte for his precious help in numerical tests concerning 
the problem with spurious fermionic modes of the WD operator and the question of 0(a) 
improvement of hadronic masses. 
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where 'D^{ip,ilj,U) is the functional integration measure, <S'ym(^) is any 
sensible lattice regularization of the Yang-Mills action and 

X ^ 

- , _ 4r 

+ip{x + ii)Ul{x){r + j^,)ipix)] + Tp{x){Mo + —)^p{x] 



E V^(^) [l E + ^m) - 4 E V^i^M + ^o] ^P{x) (2.2) 



is the standard Wilson fermion action. In the last line of eq. 1)2.2(1 we have 
introduced the usual definitions of forward and backward lattice covariant 
derivatives 

V^V(a^) = - [Uij.{x)ip{x + afi) - V'(x)] 
a 

V;%b{x) = i [^^(x) - Ulix - afi)^{x - aft)] . (2.3) 

In all our formulae the sum over flavour indices is understood. In this paper 
we restrict to the case in which Mq is taken to be proportional to the unit 
matrix in flavour space. It is, however, possible to extend the following 
arguments to encompass the case of non-degenerate quark masses. This will 
be done in refs. fM?E\. 

The piece of the action proportional to r is to the so-called Wilson 
term. It is an "irrelevant" operator of mass dimension 5, needed to solve 
the fermion doubling problem. Its presence explicitly breaks the SU{2)l 
SU{2)ji chiral symmetry of the (naive discretization of the) fermionic ac- 
tion. The symmetry is only recovered, at tree level, when Mq and a go to 
0. In higher orders of PT powers of the factor a in front of the Wilson 
term may be compensated by divergencies in loops, leading to finite or 
even divergent r-dependent contributions. This means that the formal chiral 
properties of continuum QCD are lost: the limit Mq — > does not corre- 
spond to the chiral limit and operators naively belonging to different chiral 
representations mix among themselves. A complicated renormalization pro- 
cedure is needed in order to recover perturbatively (and non-perturbatively) 
chiral symmetry in the continuum limit |^. 

These are intrinsic problems that do not depend on how fast the contin- 
uum is reached. The same kind of problems appears, in fact, in the improved 
version of LQCD |E1 13 El El , even if, in the limit a ^ 0, 0(a) corrections 
are absent from lattice estimates of physical quantities. 



6 



We are interested in computing expectation values of the type {xi / 

X2 / . . . / X„) 



{0{xi,X2,...,Xn)) ^ ^,^= (2.4) 
(r.Mo) 

= / P/x(V',^,[/)e"[^YM+^^no(xi,rE2,...,x„), 

^LQCD J 

with 0{xi,X2, • • • , Xn) a gauge invariant, multi-local and m.r. operator built 
in terms of the fundamental fields of the theory (quarks and gluons). We 
will also assume that O has in the continuum limit a Lorentz (rather, in 
Euclidean metric, an S'0(4)) covariant structure. The subscript (r, Mq) 
on the l.h.s. of eq. (|2.4|) means that the expectation value is taken with the 
specified values of the Wilson parameter, r, and bare quark mass, Mq ^. Here 
and in the rest of the paper we will always restrict to connected correlation 
functions, though this will not be explicitly indicated in our subsequent 
equations. 

The key observation of this work is that under the IZ^ transforma- 
tion (|1.3)1 the fermionic action (|2.2|) goes into itself, if at the same time 
we change sign to the Wilson term (i.e. to r) and Mq. In the spirit of the 
spurion analysis, a quick way of studying the situation is to consider r and 
Mq momentarily as fictitious fields and consider the combined transforma- 
tion 

nf = 7^5 X [r ^ -r] X [i\fo ^ -Mq] (2.5) 

as a symmetry of the lattice theory. In this situation it is easy to prove 
that any m.r. operator will be either even or odd under 7^5^. The argument 
goes as follows. Given the bare operator, Ob-, one can always assume that 
the latter is either even or odd under 7^5, because [T^s]^ = 1. Since 7^5'' 
is a spurionic symmetry of the action, the renormalization procedure will 
necessarily give rise to a m.r. operator, O, which under Ti^^ has exactly the 
same parity Ob has under TZ^. This parity will be denoted in the following 
as and called the T^-s-parity of O, while PnAO] will be referred 

to as the 75-chirality of O. 

In general O will be given by a linear combination of bare operators with 
finite or power divergent coefficients having a parametric dependence on r 
and Mq. In this linear combination operators that under TZ^ have a parity 



^Throughout the paper to hghten the notation we will drop, unless necessary for clarity, 
the dependence upon the bare coupling constant, go. Also with a little abuse of notation 
we will use the same symbol O for the operator on the Hilbert space and the function of 
the elementary fields that appear under the functional integral. 
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opposite to that of Ob (if any) will necessarily appear multiplied by mixing 
coefficients odd under [r —r] x [Mq — > —Mq]. 

This argument implies that the expectation value of a m.r. operator must 
satisfy the relation 

{0{xi,X2, . . ■,Xn)) 

where, we recall, the sign inversion of r and Mq has to be performed both 
in the action and in the mixing coefficients. This result can be proved by 
performing the change of fermionic integration variables induced by in 
the functional integral defining the l.h.s. of eq. H2.6|) . 

Prom this discussion one also concludes that the renormalization con- 
stant of O, Zq, must be even under r — > — r in order for the renormalized 
operator, O = ZqO, to have the same 75-chirality as O (and Ob)- We will 
see in sect. K-i.lHl how it is possible to extract from lattice data non-perturbative 
estimates of renormalization constants that are even in r and 0(a) improved. 

The interest of eq. (|2.6j) lies in the fact that from it one can prove that 
discretization effects associated to the presence of the Wilson term in the 
action, induce in correlators 0(a) corrections having well defined properties 
under r —r (r-parity for short). The latter are such that averages of 
Green functions computed in lattice theories with opposite values of r have 
a faster approach to the continuum limit and a better chiral behaviour than 
the two correlators separately. 

Ideas in this direction were already put forward in the literature long 
time ago in the papers of ref. ISOJI^ISl] ^■ 

To give a precise basis to the line of arguments sketched above, it is 
necessary to start by establishing the r-parity property of Mcr(r). Since 
Wilson's lattice theory (|2.H) is invariant under the spurionic transformation 
TZ^^ (see eq. ()2.5() ). the mass action counterterm a^ J2x '4'{x)Mcr{r)ilj{x) must 
be naturally chosen so as to maintain this spurionic invariance. From this 
simple but crucial remark, given the definition of 7^5, it is immediately 
concluded that the critical mass must be an odd function of r 

M,,{-r) = -M,,{r) . (2.7) 

For the reader's convenience, we illustrate how this this key property is born 
out by the definitions of Mcr that are commonly employed either in PT or 
at the non-perturbative level, where S^SB is expected to take place in the 
large volume regime. 

^We thank S. Aoki and S. Sint for having brought these papers to our attention. 



(r.Mn) 



-Mr,) 



(2.6) 
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To avoid mixing up different issues we will present the whole discussion 
in this and next section by temporarily ignoring the problems with the 
spectrum of the WD operator. We will take up this crucial question in the 
second part of this paper beginning from sect. |1J 



2.1 Critical mass in PT 

We start from the formal definition of the lattice quark propagator, S^^, in 
momentum space 

Performing the change of fermionic functional integration variables induced 
by the transformation ^\.'^^ . one gets the lattice identity 

r, Mo) = -75S(p; -r, -Mo)75 ■ (2.9) 

In the continuum inspired decomposition (where we ignore obvious logarith- 
mic factors irrelevant for this argument) 

S(p; r, Mo) = iSi(p2; r, Mo) ^iip'^r, Mo) + 0(a) , (2.10) 

the critical mass is defined as the value of Mq which solves the equation 

S2(0;r,Mo) = 0. (2.11) 

Introducing the decomposition ()2.1U() in (|2.9() . we derive the simple relations 

Si(p2;r,Mo) = Si(p2;_r,-Mo) (2.12) 
S2(p2; r, Mo) = -S2(p2; -r, -Mo) . (2.13) 

To any order in perturbation theory Si and S2 are smooth functions of Mo 
and r 7^ 0, which are defined for arbitrary values of Mq and l^l < 1. It 
follows that, if a certain value of Mo is a solution of 1)2. llj) for a given r, 
then —Mo solves the same equation with r replaced by —r. In other words, 
barring unexpected patterns of multiple solutions to eq. ()2.11|1 . we recover 
eq. (|2.7|) . in agreement with the argument given in ref. |31j . 



2.2 Critical mass beyond PT 



In the presence of S^SB M^r can be defined as the value of Mo at which the 
pion mass vanishes. This definition is implemented numerically by studying 
the large t > behaviour of the two-point correlator 



G^(t;r,Mo) 



X 



(r.Mo) 



(2.14) 
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where 

= ^TsyV' (2.15) 

is the pseudo-scalar quark density with flavour index b (b = 1,2,3). Per- 
forming in the expectation value ()2.14|) the change of functional integration 
variables ()1.3() . one flnds 

G^{t-r,Mo) = G^it;-r,-Mo). (2.16) 

As implied by reflection positivity (see sect. 12.31 below and Appendix B), we 
assume that the mass, ttt-tt, of the lightest state contributing to ()2.16|) can be 
non-perturbatively, equally well defined for positive and negative values of 
Mq and for \r\ < 1. Then to comply with the behaviour implied by eq. ()2.16|) 
in the t ^ oo limit, we must have 

m^{-r, -Mo) = m^(r, Mq) . (2.17) 

From the definition of Mcr the relation ()2.7() again follows. 

In deriving this non-perturbative result we have implicitly assumed that, 
once Qq and r are given, there always exists a unique value of Mq where 
the pion mass vanishes. Careful analyses of cutoff effects in lattice QCD 
with Wilson fermions |^ IHS] actually suggest that, depending on details 
of the lattice theory and the value of gQ, the mass of the neutral pion can 
either vanish at two "critical" values of Mq (with the charged pions being 
massless for Mq in the interval between them) ^ or never vanish and just 
have a minimum at some value of Mq (in this situation charged pions remain 
degenerate with the neutral state). In the latter case, the value of Mq at 
which m^r is minimal can be defined as Mcr(r). This definition satisfies 
eq. ()2.7() . In the former case, instead, the two "critical values" of Mq turn 
out to differ by an amount that is O(a^) in the scaling limit . In this case 
any of the two "critical" values can be chosen as Mcr(r): e.g. the larger for 
r > and the smaller for r < 0. This choice is consistent with eq. (|2.17|) 
and thus with eq. 1)2. 7|) . 

As expected, we find that for the purpose of making the subtracted Wil- 
son term in the action, i/j[—a{r/2) J2f_i ^Ji^fj. + Mcrir)]^^, a truly dimension 
five operator the quantity Mcr(r) can always be taken odd in r. 

It goes without saying that, if PCAC is used to define Mcr, an equally 
good and r-odd determination is obtained (see Appendix A), only differing 
by 0(a) terms from that based on the behaviour of the pion mass, described 
before. 

^We restrict attention to the branch of the Aoki phase closest to Mo = (see Fig. 9 of 
ref. ED). 

''We thank the referee for drawing our attention to tliis point. 
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2.3 A comment on reflection positivity 

An important question in the framework of the computational strategy we 
are proposing is whether reflection positivity j341 135j remains valid for the 
Wilson lattice theory in which the sign of r is reversed. Reflection positivity 
is a crucial property for being able to extract energies and matrix elements 
from lattice correlators, because it ensures that a consistent Hamiltonian 
formalism can be set up at finite lattice spacing (even though the explicit 
construction of a transfer matrix may be technically difficult or not avail- 
able) . 

It has been shown in refs |35| I36j that the "link" reflection positivity 
of the lattice action is a sufficient condition to define a Hilbert space of 
states with positive metric and prove the existence of a positive, self-adjoint 
operator that can be viewed as the square of the transfer matrix. We recall 
that link reflection, as opposed to "site" reflection ^37..38j . refers to reflection 
of lattice points across a time-hyperplane not containing lattice sites. 

The standard fermionic Wilson action 1)2. 2|) is known to be link reflection 
positive for all (real) values oi Qq, Mq and \r\ < 1 12^. Actually in refs. [HSl 
I36| a stronger result was proved which implies that even the general tm- 
LQCD action of eq. (|4.1|) is link reflection positive for |r| < 1. This fact 
is of special relevance for our work, where the problems with the spectrum 
of the standard WD operator are avoided by making recourse to tm-LQCD 
(see sect. HJ). 

In Appendix B for completeness we give some useful definitions and sum- 
marize a few results about reflection positivity and lattice transfer matrix. 

3 Removing 0(a) cutoff effects 

In this section we propose a way of extracting 0(a) improved estimates of 
physical quantities from lattice correlation functions, which does not require 
the knowledge of anyone of the usual improvement coefficients. 

Let us consider the expectation value of a gauge invariant, multi-local, 
m.r. operator which, as we discussed in sect.|21 without loss of generality can 
be taken to have a definite 75-chirality. In fact, while multi-local operators 
made out only of bosonic local factors are necessarily eigenstates of 7^5^, the 
general case of multi-local operators containing also baryonic local factors 
can be treated by splitting the operator into a 7^5-even and a 7?.5-odd part 
and applying separately to each of them the arguments we present in the 
following. 

We start by noting that, since Mcr(r) is odd in r, the excess quark 
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mass niq, eq. H1.2() . changes sign under Mq — > —Mq and r — > — r. The 
transformation TZ^r^ under which the Wilson action is invariant can then be 
also written as 

nf = n5x[r ^ -r] X [ruq -m,] . (3.1) 

Correspondingly, the relation ()2.6|) . which expresses the implications of the 
spurionic symmetry TZ^^ on lattice correlators, takes the form 

(O) = {-if^5lo](^o) , (3.2) 



{r,mq) 



{-r-niq) 



where, by slightly changing the notation employed in the previous section, 
we have indicated the parameters that specify the fermionic action by using, 
besides r, the excess quark mass, niq = Mq — M^r, instead of the bare mass 
Mo. 

In order to discuss the issue of 0(a) improvement we need to make 
reference to the notion of effective theory introduced by Symanzik and use 
the related Symanzik expansion of lattice correlators in terms of correlators 
of the continuum theory [3]. Schematically up to 0(a) terms, one can write 
for the lattice expectation value of a m.r. operator the expansion 

I ^ ^ cont 

l(r,mg) {mq) 
— _ ^ cont o 

+aY.{mqrr^'^^(r){0,) +0{a^). (3.3) 

The label |(^*) in the correlators appearing on the r.h.s. of eq. (|3.3|) is meant 
to recall that they are continuum Green functions renormalized at the scale 
a^^ . They are computed employing the continuum QCD action, regular- 
ized by using, e.g., a second lattice regularization with lattice spacing much 
smaller than a. Consistently, the parameter ruq on the r.h.s. is to be inter- 
preted as the continuum value of the quark mass at the subtraction point 
a^^ . Even in the limit of vanishing lattice spacing, the previous equation has 
a logarithmically divergent a-dependence. This divergency can be disposed 
of by multiplying both sides of eq. (|3.3|) by the renormalization constant of 
the operator O. The coefficients ("j ^ and r/ are finite functions of r and g^. 
All these coefficients are necessary to match the lattice correlator on the 
l.h.s. to the continuum Green functions on the r.h.s. of l|3.3)) . 

The sum over £ is extended over all operators. On ^ O, with (mass) 
dimension {d[0] is the dimension of O) 

d[Oe] = d[0] -m + l, (3.4) 



12 



having the same unbroken quantum numbers of O. The obvious positivity of 
d[0\ and d[Oe\ puts an upper bound, nmax, to the possible values of ni > 0. 
Two types of correlators will contribute to the sum over t. correlators in 
which products of local operators are inserted and correlators where O is 
inserted with space-time integrated local densities. The former are related 
to the improvement of the single local factors making up O, while the latter 
owe their origin from the need of improving the action _8_. 

A careful discussion of the structure of the Symanzik expansion will be 
given in Appendix A, where we present a detailed analysis of the improve- 
ment of flavour non-singlet axial WTI's. 

In the light of eq. (|3.2)) we want to compare the expansion (|.S.3j) with 
the similar one obtained after changing sign to r and ruq. To this end we 
observe that 

1. in the r.h.s. of eq. (|3.3() niq appears not only explicitly in some of 
the factors that multiply the continuum correlators, but, as we have 
explicitly indicated in the formulae, also through the dependence of 
the continuum action on m^; 

2. the transformation (|1.3j) is part of the chiral group and, acting on 
the continuum action, is equivalent to changing the sign of niq] since 
[7^5]^ = 1, the operators can always be taken to have well defined 
75-chiralities, P-ji^lO^]] consequently their continuum expectation val- 



ues will have a definite parity under niq 



-m, 



cont 



[Oe 



9' 
cont 



namely 



(3.5) 



cont 



Using (|3.5j) . one gets from (|3.3j) 

{0)\ . = (-l)^-«[^l[Cg(-r) - amq^Si-r)]{0) ^ ^ + 

+a ^(-l)^-5[0^1+-^(rn,)"^r?g^(-r)(0,) + 0{a') . (3.6) 

Inserting the expansions (|3.3() and (|3.6|) in eq. H3.2() and equating 0(a°) and 
0(a) terms, one obtains the relations 

C8ir)-C8i-r)=0, (3.7) 



amq[^Sir)+C^{-r)]{0) 



Oi 



{niq) 



+ 



(3.8) 



cont 
{mq) 



0. 
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We now want to show that eqs. p.7p and p.8p imply the formula (Wilson 
average - WA) 



(O) 



WA 
(rriq) 



(O) 



(r,mq) 



+ {0) 



(—r,mq). 



cont 



C^(r)(0) +0{a'). (3.9) 



We recall again than in the above equation the prescription of changing 
the sign of r must be extended to the whole r-dependence present in the 
expectation value, i.e. to the r-dependence of the action and of the mixing 
coefficients making up O in terms of bare operators. 

To prove eq. (|3.9() we first observe that from eq. (|3.3|1 , by changing only 
sign to r, one gets 



(r,mg) 



+ {0) 



(—r,mq) 



0( 



{O) 



cont 



(mq) 



+ 



(O) 



cont 



{niq) 



+ 



cont 



(3.10) 



In Appendix C we prove that simple symmetry considerations imply the 
validity of the relation 



Pn, [O] + Pn, \Oi\ +ne = l mod (2) . 



(3.11) 



Consequently all 0(a) terms on the r.h.s. of eq. (|3.10|) vanish owing to the 
relation 1)3. 8() . while the 0{aP) terms add up to the desired continuum result 
thanks to eq. (|3.7j) . thus proving our statement. 

It is worth noting that a compact way of summarizing the results ob- 
tained above is to say that the Symanzik coefficients C, are even functions of 
r, while all the others and rf) are odd. These properties immediately fol- 
low from the vanishing of the sum of terms in the last two lines of eq. ()3.1U|) 
(owing to eqs. (|3.<S|) and (|3.11|) ). thanks to the fact that in each term r and 
space-time dependences are completely factorized. 

Making use of eq. (|3.2jl in the second term of the formula (|3.9|) . we can 
rewrite the latter in the equivalent form (mass average - MA) ^ 



(O) 



MA 
(mq) 
Oi 



(O) 



{r,mq) 



+ (-1) 



(O) 



(r,—mq) 



Q'6{r){0) +0{a') 

[mq) 



(3.12) 



'We thank M. Testa for a very useful discussion on this point. 
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which says that an 0(a) improved lattice correlator can be obtained by tak- 
ing the sum with an appropriate relative sign of the two correlators computed 
within the same lattice regularization (same value of r) but with opposite 
sign of the excess quark mass, m^. 

Since, as we have seen, eqs. (|3.9p and (|3.12|) are identical, in the following 
we will often refer to either one of them as Wilson average. 

3.1 Comments 

A few observations on the general strategy we have outlined above are in 
order here. A first remark concerns the key property characterizing the 
operators to which our method is applicable. Two others have to do with 
the question of whether 0(a) discretization effects in M^r may jeopardize 
the results we have obtained and the delicate problem of taking the ruq — > 
limit in eqs. (|^ and (jXT^ . 

3.1.1 Multiplicative renormalizability and improvement 

At the beginning of sect. Ol we restricted our consideration to m.r. (gauge 
invariant) operators, O. Indeed the Symanzik expansion ^'^.'A^ . in the way 
we have written and exploited it, is only valid if O is (apart from a possible 
logarithmically divergent renormalization factor) a finite operator. If this is 
not the case, negative powers of the lattice spacing would appear in the r.h.s. 
of the expansion, together with dreadful terms arising from the interference 
of 0(a) (and higher) discretization effects with power divergencies. In this 
situation the naive counting of powers of a is not anymore valid and the 
Symanzik expansion is not immediately useful for our purposes ^. 

3.1.2 0(a) discretization effects in M^r 

While in PT, to any order in g^, M^r can be defined to an arbitrary accuracy 
in the lattice spacing, this is not so non-perturbatively. Beyond PT, Mcr can 
be determined, in the presence of SxSB, by looking at the behaviour of the 
pion masses (as discussed in sect. 12.2)) or, more in general, as the value where 
the axial currents are conserved (PCAC). In any case the critical mass can 
only be known up to discretization errors whose magnitude depends on the 
condition chosen to determine it and the details of its computational imple- 
mentation. Then an obvious question to ask is whether these ambiguities 
may invalidate the argument about improvement presented above. 

®We will take back this question in |29| . when dealing with power divergent operators, 
like the effective weak Hamiltonian, and show how one can go around this problem. 
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We show in Appendix D that, although Mcr can only be known with 
0(a) discretization errors (unless the theory is fully improved a la Symanzik, 
something which obviously we do not imagine to be doing in the present 
context), the formulae for Wilson (eq. ()3.9|) ') - or mass (eq. ()3.12|) ) - averages 
lead to 0(a) improvement as claimed above. 



3.1.3 The limit rUg ^ 

In the presence of SxSB the chiral symmetry is not realized d la Wigner 
and, as the mass goes to zero, the chiral phase of the vacuum is driven by 
the phase of the quark mass term. The same must be true on the lattice, 
thus ideally the continuum limit should be taken before letting rriq 0. As 
a result, taking the chiral limit is a numerically delicate matter (see also 
sect.EXTJ) 

In practice in order to ensure that on the lattice the chiral phase of the 
vacuum is determined by the quark mass term proportional to mq, and not 
by the Wilson term, it is necessary to work with lattice parameters satisfying 

"^g^QCD » aAqcD , (3.13) 
in the unimproved case, or possibly under the weaker condition 

"^g^QCD » a^^lcD > (3.14) 

if one is dealing with 0(a) improved correlators. When conditions like (|3.13|) 
or (|3.14|) are not fulfilled, renormalized lattice quantities may strongly dif- 
fer from their continuum limit counterpart, thus resulting of little use in 
extracting information about the underlying theory. 

To analyze the subtleties of taking the limit niq ^ in our approach 
it is convenient to first discuss the non-physical case in which there is no 
SxSB. 

• In the absence of SxSB there would be no problem in setting niq = 
and from (|3.2|1 one would get 

(O) = (-l)^^5[o](0) . (3.15) 



(r,mg=0) 



(-r,mg=0) 



This relation is telling us that expectation values of operators with Pn^ [O] = 
mod (2) are automatically 0(a) improved. Indeed all the operators Oi 
that are left over after taking the limit niq (i.e. those with ni = 0) 
in the Symanzik expansion of O are necessarily T^-s-odd (see eq. (|C.2|) with 

^°We wish to thank R. Sommer for drawing our attention to this important point. 
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d[Oe] = d[0] + 1), hence have a vanishing continuum expectation value. So 
there are actuahy no non-zero 0(a) contributions in the Symanzik expansion 
of the lattice expectation value of O if -Pr,5[0] = mod (2). On the other 
hand for operators having P-ji^ [O] = 1 mod (2) the WA yields a quantity 
which is identically zero, in agreement with the expected continuum result. 

• In the presence of SxSB the chiral phase of the vacuum matters and, 
since the theory is non-analytic at niq = 0, eq. (|H.15|) does not make sense 
anymore, rather the limit — > should be approached under the condi- 
tion or (|l-i.l4j) if correlators are 0(a) improved. In this situation a 
direct analysis of the rrig-dependence of (0)|(r,m ); based on the parities of 
the continuum expectation values under niq — > —mq (see eq. (|3.5|) ). confirms 
the validity of eq. 1)3. 12() . 

The remarks of this section suggest that in practice, for each value of 
the cutoff effects on quantities sensitive to S^SB may become quite large 
when niq goes below some threshold value. This threshold value should 
be determined by numerical experiments and will depend in general on the 
observable under consideration. As a consequence, performing a chiral fit 
of lattice data obtained at fixed gQ becomes difficult and possibly inappro- 
priate, if data at too small values of niq are included. The safe region is 
determined by the order of magnitude inequality (|3.13j) or rather, for im- 
proved quantities, (|3.14|) . It is reassuring to note that for realistic lattice 
spacings, say a < 0.1 fm, the bound (|3.14|) is satisfied down to pretty low 
values of niq, of the order of 10 Mev. 

3.2 Matrix elements and hadronic masses 

An immediate consequence of the absence of 0(a) discretization effects 
in Wilson averaged correlators is the observation that improved hadronic 
masses and matrix elements can be similarly obtained by taking averages 
of the corresponding quantities, extracted from the individual terms in the 
l.h.s. of equations like (|3.9|) (or equivalently (|3.12|) ). i.e. by averaging the 
values of the masses and matrix elements computed in separate simulations 
performed with standard Wilson actions having opposite values of r and 
fixed niq (or fixed r and opposite values of ruq). 

This remark is very important as it yields a simple and practical prescrip- 
tion to extract improved physical quantities from Monte Carlo simulations. 
It also solves an apparent paradox of the approach we propose The prob- 
lem has to do with the observation that different regularizations, as the ones 

^^We would like to thank R. Sommer for a useful discussion on this question. 
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provided by Wilson lattice actions with opposite Wilson terms, will lead at 
finite a to slightly different values for the masses of intermediate states. Now 
suppose that we consider the spatially integrated two-point Green function 
of some local, m.r. (bosonic) operator, <I>/j, with definite 75-chirality and the 
quantum numbers appropriate to create the state h {\h)) from the vacuum 
(|r2)). We have proved that to get 0(a) improved quantities, one has to 
average correlators computed with opposite values of r. In particular we 
expect to be able to extract improved masses from the average 

WA 1 r -1 

Ghhit) ^ = -\Ghh{t;r,mq) + Ghh{t;-r,mq) , (3.16) 



with 



Ghhit; r, mq) = YH^hi^^ tHi^)) , , ■ (3-17) 



Here the operator is defined in terms of the lattice reflection operator, 

0, appropriate to Euclidean metric (see Appendix B for details), so as to 
guarantee the validity of the relation (a|$Jj|/3) = (/3|$/i|a)* for generic states 
|q) and 1/3). 

If we take at finite a the limit t — 00 of Ghh{t)\^l^y we seem to run 
into a paradox. On the one hand, in fact, as i — > 00 one of the two correla- 
tors, namely that in which the lowest lying state has the lighter mass, will 
dominate over the second. On the other hand, one single correlator is not in 
general 0(a) improved by itself. Moreover at finite lattice spacing it is not 
obvious that one can extract hadronic masses from the large time behaviour 
of the Wilson averaged correlator, Ghh{t)\^^^, as the latter does not admit 
a spectral decomposition in the usual sense. 

The solution of this problem lies in the remark that for the purpose of 
getting 0(a) improvement of physical quantities via I^^'s the limits t ^ 00 
and a ^ do not really commute. In fact, if we let \t\ ^ 00 at fixed a, the 
corrections of order a\t\ are no longer small. As a result the expansion in a 
breaks down and the continuum limit cannot be taken. This simply means 
that the continuum limit a ^ must be taken first. This observation, 
however, is not immediately useful for the purpose of actually extracting 
improved masses and matrix elements from lattice data. 

To see how one should proceed in practice let us fix the attention on the 
spatial Fourier Transform (FT) 

GhhiK t- ±r, rriq) = a^Y. *)^I(0)) , • (3-18) 
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From the general principles of quantum mechanics and the existence of a 
lattice Hamiltonian for any |r| < 1 (see sect. \2.'6^ . it follows |341 135| that the 
following spectral representation (t > 0) must hold 



(3.19) 



where the £'/i_„(k; itr, mg)'s denote the energies of the (covariantly normal- 
ized) states n, k)|(±j.^m^) in the two regularizations we are comparing and 



^h,n(k;±r,mg) = (il|^>/j|/i,n,k)| 



{±r,mq) ■ 



(3.20) 



For combinations of operators and states such that i?/i^„(k; ±r, m^) is non- 
vanishing in the continuum limit, we expect for sufficiently small a 



E,,„(k; ±r, m,) = (k, m,) + , 5± = 0(a 



|i?,,4k;±r,m,)|^ = |Q(Orl^M (k,"^.)r(l + e: 
From the absence of 0(a) terms in the average 



WA 



1 r 



(rnq) 2 

and the arbitrariness of t, one gets 



GhhiK t; r, TUq) + GhhiK t; -r, 



(3.21) 
0(a). (3.22) 



(3.23) 



O(a^) 



h,n 



+ e 



h,n 



and consequently 
1 



Eh,nO^; r,mq) + [r 



i?;j°f(k,m,)+0(a^ 



(3.24) 
(3.25) 



(3.26) 



1 



|i?,,„(k;r,m,)|2 + [r--r] =\C^l^{r)f\R'i^-\k,mq)f + 0{a'). (3.27) 



The practical recipe suggested by these equations is to first extract energies 
(hence masses, £'/i^„(0; r, nig) = Mh^n{r, niq)) and matrix elements separately 
from the two terms in the r.h.s. of eq. (|3.23j) and then take the WA of the 
corresponding quantities to get rid of 0(a) cutoff effects. 

We show in Appendix E that the argument leading to eq. (|3.26|) di- 
rectly applies to the energy of baryonic states, with no need of splitting the 
interpolating operators into T^s-even and -odd parts. 
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A somewhat more elaborate argument can be made to show that it is 
possible to take the "square root" of eq. H3.22() before computing the average, 
thus getting the simpler relations 



Rh,n{^;r,m,) + [r ^ -r] = Q(r)i?^-*(k, m^) + O(a^) , (3.28) 



which say that Rh^nO^jf^^q) and Rh^nO^] —r,mq) only differ by 0(a) terms 
that are odd in r 



3.2.1 Generalization to three-point correlators 



The above line of arguments can be generalized to encompass the case of 
three-point correlators. Let us consider the double FT 



G/,B/i/(k,t,k',f';±r, ruq) 



a" ^ e"" "e 

x,x' 



ikx„— jk x 



(<l>,(x,t)i?(0)cl>t,(x',t')) 



{±r,mq) 



(3.29) 



where, like also B is a local, m.r. operator. With an obvious extension of 
the previous notations we can write for the spectral decomposition of G^Bh' 
{t >0,t' < 0) 



GhBh' (k, t, k', t'; ±r, rUg) = ^ 



-Eh,n(k;±r,mq)t E^, _^,(k' ■,±r,mg)t' 



^Eh,nO^;±r,mq) 2Eh'y{k';±r,mq 



n,n' 



•i?/i,„(k; ±r,mg)(/i,n,k|5|/i ,n ,k 



{±r,mq) 



Rh',n'i^';±r, rriq) 



(3.30) 



0(a) improvement of the correlator (|3.29|) via WA, together with eqs. ()3.26|) 
and (|3.27|) (or rather (|3.28|) ). imply the set of relations 



1 r 



{h,n,\<i\B\h' ,n' X 



{r,mq) 



+ [r 



C§{r){h,n,k\B\h',n',k' 



cont 



(rriq) 



+ 0(a2 



(3.31) 



Eqs. H3.31() prove that VFA's of matrix elements of local operators with 
definite 75-chirality between pairs of states having arbitrary three-momenta 
are 0(a) improved. 

^■^If there is only one operator that can create from the vacuum the state \h,n,]s.)\(^r,mg), 
a suitable choice for its phase is sufficient to obtain eq. 13.2811 . When there are several 
such operators, consideration of all the related two-point correlators is necessary. 
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3.3 Renormalization constants 

In agreement with the general argument about the r-dependence of renor- 
malization costants given in sect. |21 (below eq. ()3.2|) '). we exphcitly show that, 
in a mass independent renormahzation scheme, the lattice renormalization 
constant, Z^, of a m.r. local operator, B, can (and should) be defined as 
an even function of r. Furthermore we prove that Zb can be determined 
non-perturbatively with only O(a^) discretization errors. 

To define the renormalization constant of the m.r. operator, B, let us fix 
our attention on the correlator 



with <I>B a m.r. operator located at x 7^ 0, with the same unbroken quantum 
numbers and T^-s-parity as B. Schematically, mimicking the standard con- 
tinuum construction, one can get a non-perturbative estimate of the lattice 
renormalization constant of B using the formula 



where fj, plays the role of subtraction point. At the cost of introducing some 
infrared regulator, all correlators in (|3.33|) are evaluated at niq = 0. This is 
not in conflict with what we said in sect. H3.1.3|) . because the phenomenon 
of SxSB is not expected to affect in any way the value of renormalization 
constants. 

Eq. (|3.33|) automatically implies that Zs^afi^r) is even in r, because 
each correlator in (|3.33|) is 7^5-even. Any other definition of renormalization 
constant will have to be even in r, otherwise the renormalized lattice oper- 
ator, B = ZbB, could not have the same TZ^ (hence chiral) transformation 
properties as the corresponding continuum operator. 

In actual simulations one has to relax the condition niq = 0, unless 
the Schrodinger functional setup is employed. But one may define at non- 
vanishing mass a renormalization constant, ZB{afi,amq,r), still employing 
a formula like 1)3. 33() . To avoid introducing a new symbol, we gave to this 
determination of the renormalization constant the same name we used be- 
fore, in view of the fact that Zsia^, 0, r) = Zsiafj,, r). The important point 




(3.32) 




(3.33) 
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is that Zb{o,h, amq,r) is no longer even in r and near to the continuum hmit 
differs from ZB{ci^-,r) by O(amg) terms. 

An 0(a) improved r-even estimate, can be however obtained either 
from the equation 



WA 

(mq) 



WA 

[{'^B{x)^Bm. 



|x|=M-l 



{<^b{x)B{0)) 



WA 
(rriq) 



[{^B{x)<^Bm 



WA 
(mq)' 



1/2 



\X\=fl~ 



(3.34) 



or, as we did in the case of masses and matrix elements, by taking the WA of 
the unimproved estimates computed with fermionic actions having Wilson 
terms of opposite signs. In formulae 



Z\{a^^amq,r) = - ZB{afx,amq,r) + ZsiafJ-Tamq, —r) + O(a^) . (3.35) 



The arguments developed in this section apply, in particular, to the 
renormalization constants of the pseudo-scalar and scalar quark densities, 
Zp and Zgo (we recall that z^Zp^ = Z~^l [SHI El is the quark mass renormal- 
ization constant, see eq. HA.8|) of Appendix A), as well as the renormalization 
constants of the vector and axial vector currents, Zy and Zj^. 

Also the renormalization constant of the gauge coupling, Z^, is even in 
r. A simple way of seeing it is to observe that the whole r-dependence of 
Zq is through the fermionic determinant, which at mq = is even in r. 

Concretely Zg could be computed by making reference to the quark- 
antiquark potential, Vqq^ using the equations 



Zg{a/d,r)gl = g'^g{a/d,r) , 
2 r /J \ - dVqq{d/a, r) 



4tt ■ 



dd 



(3.36) 
(3.37) 



where Cp = 4/3 for the gauge group SU{3)c and Vqq can be extracted from 
the expectation value of two Polyakov loops located at spatial distance d {d~^ 
plays the role of subtraction point, fi). In a mass independent renormaliza- 
tion scheme, we can go in the chiral limit and observe that the Polyakov 
loop correlation function is even in r. The same holds for dVqq{d/a,r)/dd 
and consequently for Zg{a/d,r). At rriq 7^ an 0(a) improved r-even deter- 
mination of Zg can be obtained through an average of the type (|3.35|1 . 



22 



We close this section by remarking that, in the context of the standard 
0(a) improvement a la Symanzik, the improvement coefficients csw lEli 
and 6g (HI E] are necessarily odd under r — > —r. Similarly odd must be the 
improvement coefficients of quark bilinears IIUI l4Uj . These results follow 
from the observation that the structure of 0(a) terms in the unimproved 
Wilson theory is dictated by the spurionic invariance of the lattice Wilson 
action under 7^5''. In particular, since these terms must be invariant under 
7^5^^, they can only be canceled by similarly invariant improvement 
counterterms added to either the action or the operators. To conclude that 
Csw-, bm, bg are r-odd, it is thus sufficient to observe that they multiply 0(a) 
improvement terms in the action which are odd under TZ^ x (m^ — > —niq). 
Similarly quark bilinear improvement coefficients are all r-odd, because the 
0(a) terms they go along with have transformation properties under TZ^ x 
{niq —friq) opposite to that of the quark bilinear in which they appear. 
This argument can be extended to a generic operator with the result that 
all improvement coefficients are r-odd. 

Actually one can go one step further and make use of this kind of spu- 
rion analysis to deduce the r-parity properties of the coefficients of "wrong 
chirality" mixings. This topic will be central in the discussion of ref. j29j . 



4 Twisted mass LQCD 

In this and the following sections we want to extend the idea of Wilson 
averaging to tm-LQCD. For completeness and to set our notations, we start 
by giving a short introduction to tm-LQCD which extends the presentation 
of ref. |2ni. 



4.1 The action 

For the purpose of this paper it is convenient to begin with a slightly more 
general expression of the lattice regularized action than the ones discussed 
in refs. [I9l| I20j. We consider a fermionic action of the form 



2 
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describing an SU{2){ flavour doublet of mass degenerate quarks Besides 
r, three more real parameters, Mg,Mg and uJr, appear in eq. H4.1|l . though, 
as we shall see, only two are actually relevant. 

As for the question of reflection positivity of tm-LQCD, it was shown in 
ref. [201 that the most general form of a (link) reflection positive (bilinear) 
fermionic action is such as to include the expression ()4.1|) for all values of 
gQ, Mq, Mq and |r| < 1 (see Appendix B for more details). 

The standard Wilson action (in the notations of eq. 1)2. 2() ) is recovered 
setting UJr = 0, Mq = 0, while tm-LQCD as formulated in is obtained 
for UJr = 0,Mg = mo and Mg = —fJ^q- The spurionic symmetry of 5p~*™ 
analog of 1)2. 5(1 takes the form 

TZf ^ 7^5 X [r ^ -r] x [M^ ^ -M^] x [< ^ -M^'] . (4.2) 

It is immediately seen that under the (non-anomalous) chiral change of 
variables 

V'(x) ^ i^'ix) = T{e)-^ij{x) 

(4.3) 

tpix) ^'{x) =ip{x)T{e)-^ 
with 

T{9) = exp(i075r3/2) , (4.4) 
the fermionic action 1)4. 1|) transforms as follows 

Sy-'^[ij,ij,U] a^^V^(x)[V7M(V*+V^) + 

X fJ, 

-a- ^ V;.V^ exp[-i(w^ + ^)75T3] + Mg cos Q - Mg' sin + 

+(M^' cos e + M'q sin 0)(-i75T3)] i>{x) . (4.5) 

With an eye to the standard expression of the formal continuum Dirac ac- 
tion, one can without loss of generality use this freedom to somewhat sim- 
plify the form of the regularized lattice action (|4.H) . It goes without saying 
that in computing expectation values of the kind 1)2. 4|) . the fields in the ac- 
tion as well as those in O must be simultaneously rotated, should one be 
willing to make use of 1)4. 

^^The case of non-degenerate masses requires some special considerations and will be 
discussed in refs. 1281 1291. 
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At Mq = Mq = the action ()4.1|) is invariant under the discrete trans- 
formation [xp = (— X, t)) 



r Uo{x)^Uo{xp) Uk{x) ^ Ulixp - ak) k = 1,2,3 



V^^ : < 'lp{x) T{u)r)lQT{LOr) ^1p{xp) 
^ 1p{x) '4>{xp)T{uJrY^-fQT{uJr) 

and the continuous ones {h = 1, 2, 3, no sum over b) 

tpix) T{uir) eyip{i9bTi,/2)T{uir)~^ij{x) 



(4.6) 



'4}{x) ^ lp{x)T{LOr) ^ exp{—i0i,Ti,/2)T{LJr) 



(4.7) 



In connection with these formulae we note that Vuj^ is the "rotated" version 
of the standard parity operation {Vujr=o = "P, see eq. (|5.4|) below) and that 
the factors T and in 1)4. 7|) cancel if 6 = 3. 

The (1-point split) currents associated with the flavour transformations 
(jlT]) are 



i^{x)-^bi^T{uJr) - r]Uf,{x)'il^{x + a/1) + 



+i^{x + aA)^[7^r(w,)"^ + r]Ul{x)^{x] 



1,2 



(4.8) 



i^{x)'^[-i^,-rT{u)r) '^]U ^{x)'4}{x + ajl) + 



+i;{x + afi)j-[j^ + rT{LOrr^]Ul{x)i>{x) 



(4.9) 



Of course is a vector current, which is conserved even when the condition 



Mq = is not fulfilled. The other two currents are mixtures of 



vector and axial currents with flavour indices 1 and/or 2, whose relative 
magnitude depends on the value of the angle Ur- We will discuss in sect. 14.21 
the structure of the non-singlet WTI's of tm-LQCD. 

The action (|4.1j) is also invariant under charge conjugation 



C : <^ 



r u^{x)^u;{x) 

^ ip{x) -V'^(x)z7o72 



(4.10) 
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and the two discrete transformations 

' Uo{x) ^ Uo{xp) Uk{x) ^ Ulixp - ak) A: = 1,2,3 



V 



1,2 



(4.11) 



Charge conjugation and Vp invariances rule out action counterterms of the 
form a'^J2x'4'{x)J2^i7fj.{'^*n + ^ ii) lbT"^{x) with r any matrix in flavour 
space. In analogy with the definitions 1)4. 11() . we could also introduce the 
transformation V^. The latter, however, as soon as cj^ 7^ or Mq 7^ 0, is 
not a symmetry of the action 1)4.1(1 . 

The existence, at vanishing Mq and Mq , of the symmetries (|4.6() and (|4.7|) 
implies in particular that radiative corrections can only generate (possibly 
divergent) mass terms invariant under V^j^ and 2^^^. Consequently the crit- 
ical mass will contribute the lattice action a single term having the form 
ipMcr exp(— za;j.75r3)V' with M^ the critical mass of the standard Wilson the- 
ory. To prove this statement it is convenient to first rewrite the action (14.1(1 
in the form 

U]=a'J2 ^(^) [I E ^f^iK + + (4-12) 
+ ( ~ E ^M^M + ^cr(r)) exp(-icj,.75r3) + mq exp(iu;m75T3) 'il){x) , 

where both Wilson and mass term are twisted. To get eq. (|4.12(1 from (|4.1|1 
we have introduced besides the excess quark mass 



Mcr cos UJr = ma COS UJ„ 



also the excess twisted quark mass 



Mq — Mcr sin iOr 



(4.13) 



(4.14) 



with niq > 0. Physics is unaffected if we subject (|4.12(1 to the transforma- 
tion ()4.3|1 with 9 = —UJr (change of variables in the functional integral). In 
the new variables the fermionic action reads 



X /I 

+ ( ~ E ^/^^M + ^^cr{r)j + mq exp[i{uJm + t^r)75T"3] X{x) . 



(4.15) 



26 



We gave it the name 




because with the identifications 

mq COs(Wm + Wr) — > fUg 

mgsin(wm + Wr) IJiq 



(4.16) 
(4.17) 



eq. ()4.15|) is precisely the tm-LQCD action of ref. ^Hl- In the form ()4.15|) 
it is clear that M^r cannot depend on ujr (nor on oJm)-, because the whole 
ujr (and OJm) dependence is brought by the excess quark mass term. The 
argument also shows that the value of the critical mass is the same as for 
the standard Wilson theory. 

Going back for a while to the form (|4.12j) of the tm-LQCD action, it 
is worth noting that 1) the transformations V^^ and X^^,, b = 1,2, are 
exact symmetries if mg = for any oj^ and, if Ur = —uJm, even when 
"^g 7^ 0) provided nig is a multiple of the identity in flavour space; 2) since 
the currents given by the eqs. (|4.8|) and (|4.9)1 are conserved in the massless 
limit, they do not suffer of any renormalization. We will make use of this 
observation in sect. 16.21 

Starting from the action 1)4.12(1 . we might equally well set the phase of 
mqtpil) to zero by bringing it fully to the Wilson term with a rotation ()4.3() 
of an angle oJm- Precisely because in this basis the quark mass term has the 
canonical form, we will call it the "physical" basis. In this basis the action 
takes the form 



^'''^V'ph, V^ph, t^] = E V5ph(x) [- E 7^(V;i + V^) + (4.18) 



and, to avoid any possible confusion, the quark fields have been denoted 
by V'phiV'ph- Eq- (|4.18() shows that physics can all be extracted from lat- 
tice quantities that only depend on three real parameters, rUg, r and the 
combination 



Eq. (|4.20() shows that the fermion determinant is real and positive and 
no zero modes can occur if mq,uj ^ 0. 




uj = ujm + oJr mod (27r) . 
In particular the fermion determinant, Pi?, has the expression 

Vp = det {D^ + mg cos uj)\D^ + mg cos uj) + rn^^siv? UJ , 



(4.19) 



(4.20) 



= E [^7m(V;; + V^) - a^V;;vJ + M„(r) . 



(4.21) 
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Ending this section, we wish to remark that the tm-LQCD action (|4.18|) 
with |a;| = -7r/2 is closely related to the action first proposed by Osterwalder 
and Seller (OS) jSS] There is, however, a key difference between the OS 
proposal and tm-LQCD. While in the action of ref. 35 quarks of different 
flavours are all twisted with the same angle, in tm-LQCD quarks of differ- 
ent flavours appear combined in pairs and the twisting involves a traceless 
flavour matrix, which for degenerate quarks is conveniently chosen to be rs 
(besides, tm-LQCD is defined for any twisting angle in the 0-27r range). This 
specific pairing of quarks is crucial to guarantee the validity of two related 
properties, which are not satisfied by the OS action |421 143j : positivity of the 
fermionic determinant (see eq. 1)4.20(1 ') and independence from the twisting 
angle of all the continuum renormalized correlation functions. 

The second property can be proved in different ways. One way is to rely 
on universality arguments ^Hl and notice that for GW fermions a twist in 
the mass term of the non-singlet type, like the one we are considering here, 
can always be let to disappear from the lattice action by a non-anomalous 
chiral transformation. This is the same as saying that in the continuum 
limit physics will not be dependent on the value of the twisting angle. An 
equivalent, perhaps more direct, way is to notice that in tm-LQCD, through 
a chiral change of fermionic variables in the path integral, one can always 
bring the whole twist on the "subtracted" Wilson term (see eq. p.lHI) ). The 
latter on the basis of symmetry and power counting arguments can only mix 
with dimension 5 operators (in particular it cannot mix with iTr[F^iyFf^i^], 
because of the different transformation properties under Vp"^ (eq. H4.11() l. 
thus it appears as a truly irrelevant operator. On the contrary, if one con- 
siders a fermion regularization with a flavour singlet twisted Wilson term, 
like the one proposed in [HS], a term proportional to i J d^x Tr[F^,yF^,y](x) 
multiplied by the twisting angle will occur in the effective gauge action. In 
this situation the corresponding continuum theory will explicitely depend 
on the twisting angle, with the latter playing the role of the 0-angle of the 
vacuum j43j . 

4.2 Non-singlet Ward-Takahashi identities 

In this section we want to write down explicitly the expression of the renor- 
malized operators (currents and quark densities) entering the flavour non- 
singlet WTI's of the theory described by the general fermionic action H4.12() . 
where both the Wilson and the mass term are twisted. The results that 

^*The use of the OS formulation to solve the problem of exceptional configurations was 
first suggested in ref. |27|. 
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we present below can be derived in various equivalent ways. One possibil- 
ity is to follow the general procedure of ref. and apply it to the present 
situation. A shortcut is to start from what has been already proved in 
ref. jl9j and derive the required formulae by performing the chiral change of 
fermionic variables that brings from the standard tm-LQCD action to the 
action (gHH). 

Our formulae are written in terms of (bare) local currents and quark 
density operators. To make easy contact with the existing literature the 
relevant renormalization coefficients are expressed in terms of the (r-even) 
renormalization constants, Zy, Z^, Zp and Z50, of the local operators of 
the standard Wilson theory (see Appendix A for notations and definitions). 

The renormalized vector and axial currents can be taken to be 



y2 



a} 



Zvv{^r,^)i'lpi^llj + ZvA{^r,^)'<lJl,jil5^'4' (4.22) 

Zvv{^r,^)'ll^l^i — ^ - ZvA{^r,^)'<lJliJil5^'4' (4.23) 

Zvi^lJ^^ (4.24) 

ZAA(t^r,'^) V'T^TsyV' + -Z'Av(wr.,w)'07My'^ (4-25) 

ZAA(t^r,w) V'7M75yV' - ^AV'(Wr,w)'07Aiy'0 (4.26) 

^aV57m75yV', (4.27) 



where Zyy, ZyA, Zav, Zaa are finite functions oi Qq, r, ujr and w, given by 
the formulae 



Zvvi'-^rj ^) 

Zva{^v,^) 

ZAAi^r,^^) 
ZAvi^r,^^) 



1 



Zm{uj) 
1 

Zm(uj) 



1 



Zm('^) 

1 

Zm{lo) 



ZyZm cos uJr COS w + Za sin uj^ siu UJ 

— ZyZm sin ojr COS uj + Za COS uJr sin UJ 

ZAZm COS UJr COS OJ + Zy sin ijjy. siu ijJ 

— ZaZui sin LOr COS UJ + Zy COS LOr sin LO 



(4.28) 
(4.29) 
(4.30) 
(4.31) 



In eqs. to (|OT|) we used the definitions (|T?T^ and (|A.10|) and we set 



Zm{uj) 



zL cos^ UJ + sin^ UJ 



1/2 



(4.32) 



With reference to the currents ()4.22|) to ()4.27|) . the WTI's with the insertion 
of the renormalized (multi-local) operator 0{y) {y ^ x) take the expected 
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form (6 = 1,2,3) 



(a*ifl(x)d(y)) 



(r,mq) 



0(a) 

2rn,{P\x)6{y)) 



+ 0(a) 



with 



rhq = Zm{'-o)Z 



m, 



1 ' 



(4.33) 

(4.34) 

(4.35) 



provided we define in terms of bare quantities the renormahzed pseudo-scalar 
operators, P'', to be 



Zp V'yTsV' 



1,2 



P = Zp Zpp{u;r,^)'4'-^lb'4^ + Zpso{u:r.,Uj)-i)'4) + Zpi_{lo)\^ 



with 



Zpp{uJr-,^) 



1 



Z 



Zpi{uj) 



Zm{^) 

1 



Zrn cos uOr COS + z sin uj.,. sin a; 



ZAf(u; 
sincij 



, sin ujr cos w + cos sin oj 



Zm{uj] 



pp{amq,uj) 



(4.36) 
(4.37) 

(4.38) 
(4.39) 
(4.40) 



Equations from ()4.33|) to 1)4.40(1 need a number of quahfications. First of all, 
the set of labels {u!r,uJm,i",mq) in eqs. (|4.33() and (|4.34l) is there to remind 
that, if the general lattice action (|4.12() is used, the formulae ()4.22|) to (|4.27() 
for the currents and (|4.36j) . (|4.37j) for the pseudo-scalar quark densities have 
to be consequently employed. Secondly, the unusual mixing of ^^75^' with 
i/jip and the identity operator is due to the breaking of parity and isospin 
induced in the action H4.12() by the twisting of the Wilson and mass term. 
Finally pp{amg,uj) is a coefficient function even in uj (and odd in r) which 
admits a simple polynomial expansion in aniq. Details about the structure 
of the power divergent subtraction in eq. (|4.37|) can be found in ref. |46j . 

In conclusion in this section we have recollected for completeness the gen- 
eral formulae (equations from (|4.22|) to (|4.40() ). that can be used to smoothly 
interpolate between the situation in which the whole twist is carried by the 
Wilson term (physical quark basis with = u;, uJm = 0, eq. (|4.18|) ) and the 
tm-LQCD formulation of ref. |12|, where tOr = 0, LOm = ^ (eq. (|4.15j) ). 
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5 Improved physics from tm-LQCD 

To make the discussion of improvement transparent it is convenient to use 
the form (|1?TH|) of the tm-LQCD action. 

With the help of the formahsm developed in sect. |21it is straightforward 
to show that VF^'s of pairs of tm-LQCD correlators computed with opposite 
values of r and a common value of ruq are free of 0(a) corrections at fixed 
value of Lo. Exactly as in the case of the standard Wilson formulation, an 
identical result is obtained if correlators computed with a fixed value of r 
but mass terms of opposite sign are combined (mass average, MA) with a 
relative sign equal to the 7^5-parity of the operator O whose expectation 
value is being computed. Moreover, an analysis of the properties of the 
(unimproved) lattice correlators and derived quantities under uj — > —uj leads 
to a deeper understanding of the practical consequences of the parity and 
isospin breaking that occurs at uj ^ 0. From this analysis it turns out that 
the tm-LQCD action (|4.18|) with lv = zb7r/2 is particularly convenient for 
physical applications. 



5.1 Improving correlators: the general argument 

The proof of improvement of VF^'s (or MA's) of correlators in tm-LQCD 
exactly parallels the argument presented in sect.|5J observing that tm-LQCD 
action (|4.18)) is invariant under the spurionic transformation (|3.1|) . We only 
remark for future use that changing the sign of the Wilson term is equivalent 
to shifting the twisting angle by tt. In other words the combined transfor- 
mation [r — > — r] X [cj — > a; -|- tt] leaves the action ()4.18|1 invariant. 

The proof proceeds by showing that all the steps we went through in the 
untwisted, Wilson case remain valid also here. One starts with the Symanzik 
expansion 



cont 
(mq) 



I (r,mq) 

- — ^ cont o 

+a J2{m,rr,S,i^,r){0e) ^ ^ + 0{a') , (5.1) 

where, as in sect. 01 O is a multi-local, m.r. operator. Lattice expectation 
values in eq. ()5.1|) are characterized by r, niq and oj. The parameter uj is kept 
fixed in all the arguments of this section and the dependence on is always 
understood. As in the standard Wilson case we are concerned with the r- 
dependence of the Symanzik coefficients ? and r]. It should be noted that 
for generic values of a; more operators, O^, contribute to the expansion (|5.1|) . 
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as compared to the standard Wilson case. The new operators arise due to 
the breaking of parity and isospin induced by the chiral twist of the Wilson 
term in eq. ()4.18|) . This fact, however, does not harm the argument that 
follows. In fact, as we prove in Appendix C, the relation 

Pn, [O] + Ptz, [Oil + n, = 1 mod (2) (5.2) 

also holds in the tm-LQCD case, because the symmetry considerations upon 
which it is based remain valid for the whole set of operators Oi entering the 
expansion (|5.1j) . Thus all the formulae we proved in sect. El maintain their 
validity also here at any fixed value of In particular all the results concern- 
ing the WA improvement of correlators (eq. (|3.9j) ) and derived quantities 
(eqs. (|3.26|) . ()3.28|) and 1)3. 31() ). as well as the property of renormalization 
constants of being even in r and the formula for the 0(a) improvement of 
their determination (sect. 13. 3() remain true in the framework of tm-LQCD. 

As in the case of standard Wilson fermions, 0(a) improvement via WA 
is equivalent to the observation that the Symanzik coefficients ^ and rj 
appearing in the r.h.s. of eq. 1)5. 1|1 have definite r-parity properties 

Cg(u;,r) = Cg(u;,-r) = Cg(^ + vr,r) 

eg(u;,r) = -eg(u;,-r) = -eg(^ + vr,r) (5.3) 
Vo.i^^r) = -rj2^{uj,-r) = -rj2^{u; +7r,r) . 

In eqs. (|5.3)) the second equality of each line follows from the invariance of 
the action (|4.18jl under [r — > — r] x [lj ^ + tt]. 



5.2 The w-dependence of lattice correlators 

While working with 7^ in the action (|4.18j) solves all the problems related 
to spurious modes of the WD operator, it implies, as we repeatedly said, a 
breaking of the parity (and isospin) symmetry. It should be noted, however, 
that if the parity operation is accompanied by a change of sign of ui, the 
action remains invariant. 

It is, in fact, immediate to verify that under the combined transformation 
■p X (a; — > —uj), where (remember the definition xp = (— x, t)) 

Uo{x) ^ Uoixp) , Uk{x)^ul{xp-ak), A: = 1,2,3 



(5.4) 
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is the physical parity of the theory, the action H4.18() goes into itself. In the 
spirit of spurion analysis this invariance and the further property x (a; — > 
— Lj)]^ = 1 imply that a definite parity label, p, can be attributed to any 
multi-local operator on the basis of the relation (see also Appendix F) 

(0(^')(X1,X2,...,X„)) J"^ = 

= {-1)P {0^P){x,p,X2P,...,Xnp))['^\. (5.5) 

(r,mq) 

We remark that the renormalization of tm-LQCD is such as to respect 
eq. ()5.5() . In particular, this means that, whenever two bare local opera- 
tors with opposite transformation properties under 1)5.41) mix, the relative 
mixing coefficient must be odd in iv. Examples of this phenomenon are 
found in eqs. KTH - K^ . (Ii:^ - and (I071) . Clearly the situation 

for V here is identical to the one we had before for TZ^. From now on, for 
short, we will call (—1)^ the formal parity of the operator Notice that 

what we are calling formal parity coincides with the parity of the operator 
when parity invariance is restored (tj = and continuum limit). 

Inserting the appropriate form of Symanzik expansion ()5.1j) in the two 
sides of the relation ()5.5() . one can prove the following parity properties of 
the Symanzik coefficients under lo —oj 

Co(p) (^'^) = Co{p)(-'^'n 

^^l:]{oo,r)=^^l:l{-uJ,r) (5.6) 

In the last of the eqs. (|5.6j) (— 1)^*^ is the physical parity of the operator Oi 
in the formal continuum theory, defined by the formula (analogous to ()5.5)) ) 



cont 



{Oe{xi,X2,...,Xn)) ^ ={-lf^{Oe{xip,X2P,---,Xnp)), ■ (5.7) 



cont 



(mg) 



From eqs. ()5.fi|) one can derive a number of interesting consequences con- 
cerning the nature and the magnitude of parity violating terms in lattice 
correlators. 

First of all, we observe that at a; = (no twisting), where the physi- 
cal parity ()5.4)) is an exact symmetry of the lattice theory, the last of the 
relations 1)5.6(1 implies 

riS':\0,r) = i-ir+P^7^S';\0,r). (5.8) 
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As expected, we find that only operators, Oi, with the same parity as 

(Pi = P mod (2)) can contribute to the Symanzik expansion of the latter, 

because p + pi = 1 mod (2) implies ??o^'*'^ ''') ~ 

If to ^ 0, parity is broken and also operators, Oi, with p£ = p + 1 mod (2) 
may contribute. These terms are odd in lo and, consistently with universality 
(they vanish at a; = 0), appear in the Symanzik expansion multiplied by an 
explicit factor of a. 

The results of the whole analysis can be summarized by saying that 

1. all physically relevant quantities are contained in the w-even parts of 
lattice correlators (see the first of eqs. (|5.6|) ). 

2. w-odd parts of lattice correlators are 0(a), 

3. contributions from continuum operators with p£ = p+1 mod (2) enter 
lattice correlators only at 0(a) or higher. 

5.3 The Lij-parity of energies and matrix elements 

Exploiting the invariance of the tm-LQCD action 1)4. 18() under the V x (oj ^ 
—to), it is possible to study the transformation properties under uj — > —a; 
of energy eigenvalues and matrix elements of m.r. local operators between 
pairs of eigenstates of the transfer matrix. Prom this analysis it follows that 
hadron masses and matrix elements between zero momentum eigenstates 
have definite parity under uj — > —uj. 

Furthermore we show that one can introduce a parity label for the eigen- 
states of the transfer matrix that turns out to coincide with the usual parity 
quantum number when parity is restored as a symmetry of the theory, i.e. 
either at w = or in the continuum limit. This notion is interesting in itself 
as it allows to control the magnitude of parity violating contributions in 
lattice correlators and proves to be useful in discussing 0(a) improvement 
in the special case uj = ib7r/2. 

5.3.1 Matrix elements, energy eigenvalues and masses 

It is not difficult to analyse the transformation properties of energy eigenval- 
ues and matrix elements under uj — > —uj and prove that masses are w-even 
quantities. 

Let us start by considering the two-point correlation function 
Ghh{^,t;uj,r,mg) = ($h(x,t) ^>j,(0,0)) 



(5.9) 
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where is a m.r. local operator with definite formal parity, (— The 
spatial FT of dEH) 



GhhiK t; u;, r, rrig) = ^ e^''-"($ft(x, t) $[(0, 0)) 



(^) 

[r,mq) 



(5.10) 



admits the spectral representation [t > 0) 

g--Bh,„(k;w,r,m,)t 



\Rh,nO^;uj,r, niq) 



where 



and n, k) 



iih,„(k;w,r,mg) = (17|$/,(0, 0)|/i, n, k) 



{r,mq) 



(5.11) 



(5.12) 



l(r,m,) 



is the eigenstate of the transfer matrix, T{uj,r,mq), that 



belongs to the eigenvalue £'/i^„(k; a;, r, mg), according to the relation 



15 



T{u;,r, mq)\h,n,'k) 



(r,mg) 

The key formula is the relation 

GhhiK t; id, r, niq) = Ghh{-K U -w, r, rUg 



{r,mq) 



(5.13) 



(5.14) 



which follows by performing in the functional integral defining the expec- 
tation value in (|5.1flj) the change of integration variables induced by the 
parity transformation H5.4|) . Eq. 1)5. 14() implies similar relations for energy 
eigenvalues and the modulus of matrix elements, namely 

Eh,n{^; io, r, mq) = £'fe,„(-k; -w, r, vriq) . (5.15) 
\Rh,n{^;^,r,mq)\ = |i?ft^„(-k;-u;,r,mg)| , (5.16) 

At k = eq. (l5T5|l becomes the announced result about the invariance of 
masses under a sign change of uj, i.e. 



M/i,„(a;, r, ruq) = Mh,n{-^, r, niq) . 



(5.17) 



In Appendix E we sketch how this result can be extended to the baryonic 
sector of the theory. 



Actually a transfer matrix has been constructed only for \r\ — 1. For |r| < 1 one 
should make reference to the transfer matrix over two lattice spacings (see Appendix B). 
For simplicity we shall ignore this subtlety in the following. 
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5. 3. 2 The parity of the transfer matrix eigenstates 



As we already remarked, the symmetry under parity of the (continuum 
and) standard Wilson theory is replaced in tm-LQCD by the invariance of 
the action under the transformation V^^ = Vx(uj^ —uj). 

Since at w 7^ parity is broken, a natural question to ask is how the 
parity operator, V, acts on the eigenstates of the transfer matrix. The 
answer is surprisingly simple and it is contained in the formula 



V\h,n,k) 



, =r]h,n\h,n,-k) Vh,n = '^- (5-18) 

{r,mq) (r,mq) 



This equation, which is proved in a costructive way in Appendix F, allows 
us to uniquely associate a sign, r]h^n = il, to each eigenstate of the transfer 
matrix. With a little abuse of language we will simply call it the "parity" of 
the state, in view of the fact that at a; = the above formula precisely yields 
the usual definition of parity. Eq. I|5.18|) implies for the matrix elements of 
a m.r. operator B with definite formal parity (— 1)^^ the important relation 



{h,n,k\B\h',n',k') =r]f^ f^,^,{h,n,-k\B\h',n',-k') (5.19) 



(r,m,) 



{r,mq) 



where 

rihn,h'n' = Vh,n{-'^T''Vh',n' ■ (5.20) 

Eq. ()5.19|) confirms that the label we have introduced coincides in value with 
the usual parity quantum number of the continuum states. In fact, at zero 
spatial momenta, k = k' = 0, the matrix element (/i, n, 0|-B|/i', n', 0)||^^ ^ 
is odd in uj if the product r]h,n'nh',n' does not match the formal parity of B, 
thus it vanishes at a; = and by universality it must also do so as a —> for 
generic lo. This means that for sufficiently small a it is 0(a) or higher. If 
instead rj^^ = +1, eq. 1)5. 19(1 at vanishing three-momenta simply tells us 
that the matrix elements of B are even in lo, thus not necessarily vanishing 
in the continuum limit. 



5.3.3 Deciding the parity of transfer matrix eigenstates 

We want to schematically describe how one can in practice decide what is the 
parity label, as defined by eq. 1)5. 18() . of a given transfer matrix eigenstate. 

Let and <I>2 be two local, m.r. operators with the same unbroken 
lattice quantum numbers and definite opposite formal parities, r/i = (—1)^^ 
and r/2 = (—1)^^ (??i = — ??2) under V. As a concrete example one may 
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think of the currents in eqs. ()4.23|) and ()4.25p . Consider the four correlators 
(i >0) 



(^) 

{r,m,) 



i,i = 1,2. 



(5.21) 



Performing the change of fermionic integration variables induced by the 
transformation ()5.4|) . it is easily proved that the correlators with i = j are 
even in co, hence of 0(1) in a, while the correlators with i ^ j are odd and 
by universality of 0(a). The spectral decomposition of the correlators Gij 
reads 



Gij (t; w, r, TUq 



2Mft„(a;,r, nig 



(^) 

{r,mq) 



, (5.22) 



where the states \h,n) 



16 

l(r,mg) 



are characterized by the set of conserved 



quantum numbers h and the integer n which labels them in order of increas- 
ing mass. 

Consider now the state lying at level n and compare the quantities 



{r,mq) 



and 



{h,n\<P2m 



(^) 

(r,mg) 



(5.23) 



that appear in the spectral expansions (|5.22|) with i = j. Since the cross 
product (il|<l>i|/i, n|$2l^)l(^m ) ^® necessarily of 0(a), as this quantity 
(or its c.c.) is the n*'' coefficient in the expansions ()5.22|) with i ^ j, we find 
that for each eigenstate one must have either 



{h,n\^i\n)\ 



{r,mq) 



or 



{h,n\<^2m 



{r,mg) 



> 



> 



(r,m.q) 



{h,n\<i>i\n)\ 



{r,mq) 



(5.24) 



(5.25) 



The reason is simply that in order for the four correlators 1)5. 21() to have the 
right "order of magnitude" in a, one has to have either 



or 



{h,n\<^i\n) 



{h,n\^i\n) 



(r,mq) 



{r,mq) 



0(1) and {h,n\^2\^) 



0(a) and {h,n\^2\^) 



{r,mq) 



{r,mq) 



0(a) 



(5.26) 



0(1). (5.27) 



'We drop the momentum label k = as it is of no relevance for the present argument. 



37 



In case the inequality (|5.24|1 holds, with r]h^n the parity of the state \h,n), 
we will set 

Vh,n = i-ir ■ (5.28) 
If viceversa the inequality (|5.25|) is satisfied, we will set 

Vh,n = (-ir • (5.29) 

Notice that the conditions (|5.24j) and (|5.25j) are not difficult to analyze nu- 
merically in the scaling region, as one is comparing 0(1) with O(a^) quanti- 
ties. For the same reason, sufficiently close to the continuum limit the value 
assigned to rjh^n will be independent from the choice of the operators $i and 
$2 and the details of their lattice discretization. This must be so, if T]h^n has 

to be interpreted as the parity label for the state n)|^^ y 

We end with an observation about isospin. Although tm-LQCD breaks 
isospin, "partial" isospin labels that become conserved quantum numbers in 
the continuum limit can be attached to transfer matrix eigenstates. This 

1 2 

possibility results from the fact that the transformations Vp can be written 
in the form 

= p X , (5.30) 

where 

: } (5.31) 

i V^ph^ ^phe-^-^1-2/2 

are discrete isospin transformations We remark that the combined trans- 
formation T^'^ X (w — > —io) sends the action (|4.18|) into itself, just like the 
product V X (uj ^ —to) does. 

Eq. (|5.30|) proves our statement. In fact, as soon as a parity has been 
attributed to a state, definite values for the isospin labels associated to 
the discrete transformations T^'^ are in turn uniquely determined, as the 
transformations Vp'^ are exact symmetries of the action (|4.18|) . 

5.4 Some observations about the practical use of tm-LQCD 

In this section we want to address two important issues concerning the 
practical feasibility of numerical simulations within a lattice regularization 
which, like tm-LQCD, besides chiral symmetry, also breaks parity. The first 
has to do with the problem of taking the limit niq — > and the second 

^^While the transformations Vp^ have the form (14.1111 in all the quark bases we have 
considered, V and T^'^ have only been reported in the text in the physical basis. 
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with the question of the scahng behaviour and the numerical importance of 
contributions with "wrong parity" to lattice correlators. 



5.4-1 The limit nig 

The subtleties of the limit nig — > we have emphasized in the case of 
Wilson fermions show up also in the case of tm-LQCD. Their discussion 
does not reveal any new feature here and we refer the reader to sect. 13.1.31 
In particular (in the presence of SxSB) to ensure that the chiral phase 
of the vacuum is determined by the quark mass term proportional to rUg, 
and not by the twisted Wilson term, it is necessary to work with lattice 
parameters satisfying the bound (|3.13l) in the unimproved case, or possibly 
under the weaker condition (|3.14j) if one is dealing with 0(a) improved 
correlators. When these order of magnitude conditions are violated, cutoff 
effects on lattice observables may become numerically large. As discussed in 
sect. 13.1.31 it remains the fact that the continuum limit should in principle 
be taken before the chiral limit. 



5.4-S "Wrong parity" contributions 

The discussion of sect. 15.31 is relevant for the analysis of actual simulation 
data, where energy eigenvalues and operator matrix elements are evaluated 
(at finite lattice spacing a, usually for nig ^ 0) by parameterizing lattice 
expectation values in terms of a spectral representation. As a consequence of 
the breaking of parity at a; 7^ 0, contributions from states with both values 
of the parity label we have introduced above will appear in the spectral 
representation. 

To illustrate the situation in a simple instance and show how one can 
possibly control these effects, let us consider the typical case of a two-point 
correlator with operators inserted at zero three-momentum 

G$H(t; cu, r, rrig) = ^(^.(x, t)E\o, 0)) J"^ (5.32) 

with $ and H m.r. local operators having equal formal parity. Prom all we 
said in sect. 15.31 it should be clear that G<i>h is even in lo. As for the matrix 
elements of the local operators between the vacuum and the eigenstates of 
the transfer matrix, there are two possibilities: 1) if the intermediate state 
in consideration has the same parity as the operators ^ and H, both matrix 
elements are uj-even and 0(1) in the continuum limit; 2) if the intermediate 
state has the opposite parity, the two matrix elements are u-odd and vanish 



39 



in the continuum limit at least as a. Thus "wrong parity" contributions, 
necessarily appearing in pairs, lead to 0(a-^) (or higher) effects. 

Arguments of this kind can be constructed for more complicated w-even 
correlators, as for instance the three-point correlators with two of the op- 
erators inserted at zero three-momentum. As intermediate state projectors 
occur in pairs in the spectral representation, contributions from states of 
"wrong parity" (compared to the naively expected one), necessarily con- 
tain an even number of w-odd matrix elements. Thus again states with 
"wrong parity" lead to contributions of O(a^) or higher. When operators 
are inserted at non-vanishing three-momentum, the situation is governed 
by eq. (|5.19)) . No new special problems arise due to the breaking of parity. 
Though a number of specific results could be derived on the subject, we will 
not pursue any further their discussion in this paper. 

Notwithstanding all the control we have on "wrong parity" contributions, 
one should not forget that, if intermediate states with naively unexpected 
parity come into play, they cannot be ignored altogether, because some of 
them may have masses smaller, or not much larger, than the state(s) of 
interest for the problem under study. Thus, for some range of values of the 
relevant time separations, the contributions from states with "wrong parity" 
can be numerically non-negligible. In this situation, if one does not want 
to end up with intolerably large lattice discretization errors on the final 
result a disentangling procedure of the type one is used to in the case of 
mixings among low- lying and excited states has to be carried out. This task 
requires the evaluation of several correlators differing among each other by 
the formal parity of the interpolating operators. 

6 A special case: cu = ±7r/2 

The choice to = ib7r/2 (full twist) in the tm-LQCD action 1)4. 18() is especially 
worth mentioning, because all quantities that are even under oj — > — w are 
0(a) improved with no need of any averaging. This result is a consequence 
of the fact that for the particular value cj = 7r/2 (or uj = —tt/2) of the twist- 
ing angle a sign inversion of the latter is equivalent (mod (27r)) to a shift by 
TT. This operation is in turn the same as inverting the sign of r. The quan- 
tities that are even in uj are hence also even in r. In this situation the two 
terms entering the Wilson average are just identical, therefore performing 
the average is unnecessary to get up 0(a) improved lattice quantities. 

Non-large 0{a^) errors on physical quantities extracted from a;-even correlators are 
however perfectly tolerable. 
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The list of physical quantities that are automatically 0(a) improved in- 
cludes among others the following important physical quantities (see sect . 15. 3)) 

1. hadronic masses 

2. matrix elements {h,n,0\B\h' ,n' ,0)\^^'^'^^ , of operators whose formal 
parity is equal to the product of the parities of the external states. 

Actually one can prove more than that. One can prove that taking appro- 
priate linear combinations of quantities computed with opposite values of 
the external momenta leads to 0(a) improved lattice data. This observation 
is very important as it shows that 0(a) improved energies and matrix ele- 
ments can be obtained from simulations carried out within a single lattice 
regularization, i.e. without the need to consider averaging over simulations 
performed with opposite values of r. 



6.1 0(a) improvement of energy eigenvalues and matrix elements 



The 0(a) improvement of averages of energy eigenvalues with opposite mo- 
menta at, say, uj = tt/2 follows from the chain of equalities 



, r, nin 



Eh,nO^'-, 



vr 
2' 



-r, nin 



h,n\ 



-k; 



TT 

2' 



-r, rria 



(6.1) 



where the first equality is a consequence of the fact that (co ^ lo + tt) x (r ^ 
—r) = 1 and the second of eq. (|5.15() . The same result could be obtained 
by using the spurionic symmetry (only valid at uj = ±vr/2) P x (r — > — r). 

At this point it is enough to observe that the WA formula for energy 
eigenvalues (analogous to the eq. (|3.26j) valid in the Wilson case) 



1 
2 

becomes 

1 
2 



+ [r 



] =Efr-\k,m,) + 0{a' 



TT 



Eh,nO^; -,r,mq) + [k 



-k] 



Efr^^k,m,) + 0{a' 



(6.2) 



(6.3) 



Concerning the matrix elements of local operators between pairs of trans- 
fer matrix eigenstates at non-zero spatial momenta, we recall that we only 
have the relation (|5.19jl from which automatic 0(a) improvement does not 
follow. Although improvement can always be achieved via WA, at the par- 
ticular value of the twisting angle we are considering in this section there 
is the possibility of getting 0(a) improved lattice matrix elements without 
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the need of carrying out two separate computations. Consider, in fact, the 
hnear combination 



{h,n,k\B\h',n',k') 



(a;) sym 
(r,mg) 



(6.4) 



1 



{h,n,'k\B\h',n','k' 



{r,mq) 



+ r]hn,h'n'{h, n, -k\B\h',n', -k') 



(^) 

{r,mq) 



Thanks to eq. (|5.19)) . this particular quantity is symmetric in uj, thus 0(a) 
improved at uj = ib7r/2. From the extension of eq. (|3.31jl to the tm-LQCD 
case with, say, uj = 7r/2 one can make contact with the continuum, obtaining 

(f )sym 



{h,n,k\B\h',n',k') 
1 



(r,mq) 



(6.5) 



{h,n,k\B\h',n',k') 



(f) 

{r,mq) 



+ (/i,n,k|5|/i',n',k'; 



{-r,mq) 



C^{l,r){h,nMB\h',n',k!) 



cont 



+ 0{a' 



where in the first equahty the spurionic symmetry x (r — > — r) has been 
used again. 

Symanzik expansions too simphfy noticeably at w = ib7r/2. Combining 
eqs. (|5.3|) with (|5.6|) . one obtains for an operator of formal parity (—1)^ 
(eq. 1)5. 5() ) the set of relations 

Co(rt(V2,^) = Co{P)(-V2,r) = Co(rt(V2,^) 



eg;:;(./2,r) eg:;(-vr/2,r) -eg:;(V2,r) 



(6.6) 



<(./2,r) ^^-P (-!)-+-< (-vr/2,r) 



^^•i^-(-ir-<(vr/2,r). 



From the above equations one concludes that the coefficients ^^(pj (±7r/2) 
are vanishing and that the only operators, Oi, which can contribute to the 
Symanzik expansion of (0*^^^)|(±^/2,m ) have = p+1 mod (2), i.e. are those 
which are multiplied by Symanzik coefficients odd under a sign inversion of 
the twisting angle. This analysis confirms that the whole mass spectrum and 
all the matrix elements of local operators that are even under sign inversion 
of the twisting angle are automatically 0(a) improved. 

The same results could be obtained by noticing that the action (|4.1<S|) 
with UJ = ib7r/2 is invariant under the combined transformation 



7^5 X P X (mg 



-nia 



(6.7) 



""^^This transformation is the product of TV^ (which is valid for any ll>) times Vx{r ^ — r) 
(which only holds at cu — ±7r/2). 
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Analogously to what we did in sects. 151 and 15. 11 to get eqs. (|6.6() it is enough 
to equate the Symanzik expansions of the two sides of the equation 



(_l)P^5[0]+P(0(f)(xpi,Xp2,...,Xp„)) 



2 ' 

{r,mq) 

(§) 



{r,-mq) 



which is obtained performing in the functional integral the change of in- 
tegration variables induced by the transformation TZ^ x V. At this point 
one has to make use of the parity and chirality properties of the continuum 
theory as well as the relation 

We wish to end this section by explicitly noting that, remarkably, the 
critical WD operator appearing in the fcrmion action (|4.18j) at liJ = ib7r/2 

D^tr^^'^ = \T.^,{K + ± ^75T3 {a'- Y: K"^, - M.r{r)) (6.9) 

is anti-Hermitian, so its spectrum is purely imaginary. Consequently the full 
WD operator cannot have any vanishing eigenvalue as soon as ruq ^ (see 
eq. (lOnj) at a; = ±7r/2). 

6.2 An application: the computation of 

As an example of the peculiarity of the choice uj = ±7r/2, we wish to discuss 
the computation of the pion decay constant, Ft^. It was already pointed out 
in refs. |221ll9j that within the fully twisted formulation of LQCD can be 
obtained from a conserved lattice current (i.e. one with trivial renormaliza- 
tion properties) and without 0(a) cutoff effects, provided the lattice action 
is improved a la Symanzik. Here we restate this result in our notation and 
show that, in fact, in tm-LQCD with uj = ±7r/2, there is no need to improve 
the action. 

In the formal continuum theory is defined by the equation "^^ 

^ -, cont ^ 

{n\Al\TT) ^ =m^V2F^, (6.10) 



^"strictly speaking eq. 16.101 is not written in a completely correct way. To comply 
with eq. 13.51 (and always taking > 0) a factor equal to the sign of rUq should appear 
in the r.h.s. As it is, eq. I|6.1()^ is valid for niq > 0, given our conventions concerning the 
relative sign of the kinetic and mass terms in the fermionic action. 
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with Aq the temporal component of the renormahzed axial current, and can 
be extracted by evaluating at large times, e.g., the ratio (t > 0) 



2 a3Ex(i^(x,t)pi(0)) 



(6.11) 



where is the pseudo-scalar quark density operator. In H6.11() we did not 
multiply by its renormalization constant, because all Zp factors would 
anyway drop out in the ratio. The isospin index has been set equal to 1, 
which corresponds to a charged pion. The reason for this choice is that in this 
way other contributing intermediate states are substantially heavier than 
the charged pion, including states with "wrong parity" (see the concluding 
remark in sect. 15. 4T^ . 

Among the possible choices for the lattice (discretization of the) axial 
current to be used in eq. (|H.11|) . the most naive one is to take the current 



Al of eq. (|4.25|) . which at cOr 



UJ 



-it/2 simply becomes 



Zv'fpph (x) —7/^75 V'ph (x) 



(6.12) 



The operator has a definite (positive) T^-s-parity (in the sense of the 
definitions H2.6|) or H3.2|) '). At the cost of a preliminary determination of Zy, 
an 0(a) improved estimate of F^^ is then obtained with no need of averaging. 
We are, in fact, dealing with zero three-momentum correlators, which are 
even under a sign inversion of the twisting angle (sect. [^H.2|) ). 

A remarkable fact about the choice ib7r/2 is that it is possible to get an 
0(a) improved estimate of Ft^ which requires neither Wilson averaging nor 
the computation of any renormalization factor. Let us, in fact, consider the 
1-point split current [22l 



Alix) 



i-pt 



1 



'^ph{x)—jfj.'y5U^{x)iljpi,{x + afi) + 



+ipph{x + a/i)y7^75C/^(x)V'ph(2;) 
'^ph{x) — U^{x)il)ph{x + afi) - 



+ 



i^phix + afL)^Ul{x)'il^ph{x) . (6.13) 



The current ()6.13|) has similarly to the current (|6.12() a definite (positive) 
7^5-parity. It is obtained from the eq. ()4.8() with 6 = 2 by setting cUr = tt/2. 
With this choice the transformation ()4.7|) is an axial transformation with 
isospin index equal to 1. Since, as we remarked in sect. 14.11 this current is 
exactly conserved at = 0, its renormalization constant is equal to 1. 
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The net result is that with the choice lo = ib7r/2 of the twisting angle in 
the action ()4.18|) and the use of the (conserved) 1-point split current (|6.13|) 
one can get an automatically 0(a) improved determination of F^^ which 
furthermore does not require the computation of any renormalization con- 
stant. This is certainly a remarkable result, especially if one observes that 
no problems can occur with spurious fermionic modes, while the required 
computational effort is minimal. 

7 Conclusions and outlook 

In this paper we have shown that in LQCD with Wilson fermions it is possi- 
ble to improve the approach to the continuum limit and the chiral behaviour 
of correlation functions of gauge invariant m.r. operators by taking arith- 
metic averages of pairs of correlators computed in theories regularized with 
Wilson terms of opposite sign and identical values of niq (Wilson average), or 
equivalently the appropriate linear combination of the correlators computed 
in a given regularization but with opposite values of ruq (mass average). 
Improved hadronic masses and matrix elements can be similarly obtained 
by taking the Wilson average of the corresponding quantities computed in 
the two different regularizations, or the appropriate linear combination of 
quantities computed with opposite values of niq. 

To avoid the problems related to the nature of the spectrum of the WD 
operator we advocate the use of tm-LQCD for the actual computation of 
the correlators taking part to the averages. We have shown that all the nice 
cancellations of 0(a) terms that we find in the standard Wilson case extend 
to tm-LQCD with mass degenerate quark doublets. Moreover, it turns out 
that, if the tm-LQCD formulation with |a;| = 7r/2 is employed, all hadronic 
energies and operator matrix elements can be evaluated with no 0(a) lattice 
artifacts and without making recourse to W^^l's, but taking (appropriate) 
linear combinations of quantities evaluated with opposite external momenta. 

A striking, though still partial, confirmation of the validity of the ap- 
proach outlined above and of the remarkable properties of tm-LQCD at 
|a;| = 7r/2 comes from the recent work of ref. |44j . where a (quenched) study 
of the scaling behaviour of the pseudoscalar meson decay constant and the 
vector meson mass was carried out down to small values of (3. Lattice results 
for these quantities at fixed pseudoscalar meson mass (~ 700 MeV) show 
surprisingly small cutoff effects from (3 = 6.2 to /? = 5.85. Furthermore a 
nice analytic check of automatic 0(a) improvement of the pseudoscalar me- 
son mass and decay constant in tm-LQCD at = 7r/2 has been obtained 
in lattice chiral perturbation theory by the authors of ref. |45j . 
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Our presentation of tm-LQCD, which follows a perspective complemen- 
tary to that of ref. provides a simple understanding of the many sim- 
plifications that have been previously observed to occur at = 7r/2 [T^ 
I2U1 1211 1^ . in the renormalization and 0(a) improvement of the matrix el- 
ements of certain local operators, such as the isotriplet axial current or the 
isosinglet scalar density. 

All the results we have obtained are based on the observation that Mcr(r) 
is odd under r — r, the property (|C.8|) enjoyed by lattice and continuum 
correlators (see Appendix C), as well as the behaviour of 0(a) terms under 
the spurionic symmetry TZ'^^ (eq. (|3.H) '). Considerations of this kind can 
be generalized without any essential modification to the case of mass non- 
degenerate quark doublets [IHl- The point will be further discussed in a 
forthcoming paper [22], where we deal with the strategy for computing the 
matrix elements of the (CP-conserving) effective weak Hamiltonian. 

Before briefly illustrating the potential of the method we have presented 
in this paper for further physical applications, we wish to comment on two 
important, related aspects of our approach, namely its computational cost 
compared to other schemes and the problem of how to most effectively per- 
form the statistical analysis of available simulation data. 

7.1 Computational cost 

The computational cost of the approach we propose is in the worst case (i.e. if 
averaging is required) twice the typical one for simulations with unimproved 
standard Wilson fermions or tm-LQCD pUl UHl HTlllHl HH] . 

The gain offered by our strategy is, however, quite significant in practical 
terms, because the systematic removal of 0(a) lattice artifacts, allowing a 
better control of cutoff effects, gives the possibility of working with coarser 
lattices than in the absence of improvement. This may lead to a dramatic 
reduction of the computational cost especially in the unquenched case, where 
CPU times have been estimated to increase like with /c ~ 7 In tbe 
quenched case, assuming "ideal" critical slowing down, one finds a similar 
behaviour with k = 5. It should be also mentioned that the efficiency 
of algorithms for simulations with dynamical fermions of the multi-boson 
type jS3 is expected to significantly benefit from the absence of the SW 
term (|1.1() in the lattice action. 
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7.2 Statistical analysis of simulation data 

A little thinking reveals that the most useful way to make use of simulation 
data depends on whether the fermionic determinant (|4.2U|) (and hence the 
partition function) is invariant under r —r (or equivalently under the 
shift a; ^ u; + vr) or not. 

If the determinant is not invariant under r — > — r, the computation of 
Wilson averages requires to carry out two independent simulations lead- 
ing to two different Markov chains of gauge configurations. The statistical 
analysis of the corresponding two sets of estimators, {Xj^r} and {Xj^^r} 
(j = 1, . . . , A'meas), of the physical quantity X must be carried out sepa- 
rately. Only after this step has been completed, the Wilson average of the 
two mean values and the associated statistical error can be computed. 

When the determinant is invariant, the whole functional measure is in- 
variant too and only one Markov chain of configurations needs to be gen- 
erated. In this situation the optimal procedure is to perform the Wilson 
average at the level of the individual (possibly "jackknived") statistical es- 
timators, Xj^r- This is done by first computing the averages 

Xf ^ = i {Xj^r + Xj,_r) , j = !,■■■, iVmcas , (7.1) 

and then going on with the standard statistical analysis of the N^^as Wilson 
averaged quantities X^^. 

As we already noticed, in all cases the extra computational effort neces- 
sary when averaging is required is not larger than a factor ~ 2 as compared 
to the case of simulations with unimproved fermions. Most interestingly, 
this factor reduces essentially to 1 in the case of simulations with \lo\ = 7r/2, 
because Wilson averaging is unnecessary. 

From this discussion it emerges that the best strategy to deal with the 
unquenched theory is to "decouple" sea and valence quarks, by having all the 
sea quarks assembled in flavour doublets and regularized with twisting angles 
equal to ±7r/2, while for valence quarks each flavour is introduced with the 
chiral phase of the Wilson term chosen so as to minimize the complexity 
of the mixing pattern of (valence) quark composite operators. Unlike what 
is usually done, here we are proposing to regularize sea and valence quarks 
with Wilson terms of not necessarily identical structure [20] • 

7.3 Outlook 

The last considerations bring about the question of the computation of the 
matrix elements of the various (valence) quark operators that are of interest 
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for physical applications, like, for instance, the chiral order parameter and 
the matrix elements of the effective weak Hamiltonian, TCcS- In particular 
for the CP-conserving sector of Ti.cs we will show in j29j that, with a clever 
choice of the twisting angle entering the Wilson term of the different valence 
quarks, it is possible at the same time to have 0(a) improvement and to get 
rid of all the "wrong chirality" mixings. The virtue of tm-LQCD in simplify- 
ing the mixing pattern of CP-conserving four fermion operators in both the 
AS' = 2 and the AS = 1 sector was already observed in refs. ^lOj and pij . 
respectively. Besides achieving 0(a) improvement, we will prove j2Hl 1201 
that mass non-degenerate sea quarks can be regularized so as to also have 
the resulting fermionic determinant real and positive. 

We would like to close this list of applications by mentioning another 
important case, where the strategy advocated in this work can be of help. 
This is the computation of the decay constants, semi-leptonic form factors 
and masses of heavy-light mesons, which has been recently shown to be 
feasible with relativistic 0(a) improved c and, even, b jSSj quark, by 
making recourse in the last case to a clever use of finite size scaling techniques 
to bridge the gap between the largely different scales present in the problem. 
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Appendix A - Improving axial Ward— Takahashi identities 

In this Appendix for illustrative purposes we discuss in full detail the im- 
provement of the flavour non-singlet axial WTI's in the framework of the 
standard Wilson formulation. 

From the action 1)2.2(1 one derives the lattice identity 

{dlA^Jx)P^{0)) = (A.l) 

= 2Mo{P\x)P\0)) + {X\x)P''{0)) 



(r,Afo) 



(r,Mo) 



where d* is the discretized backward derivative The lattice expressions 
■^^For notational simplicity we write formulae with derivatives inside expectation values. 
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of d*A'^^ and X^, whose exphcit form we do not need here, can be found in 
ref. H] . A^^ and are the axial current and the pseudo-scalar quark density 
(eq. H2.15() ). respectively. To be concrete the inserted operator in HA.1|) was 
taken to be the pseudo-scalar quark density. It would not be difficult to 
generalize the arguments that follow by considering the insertion of a more 
general m.r. operator with definite 75-chirality. 

In ref. |4| it is shown how one can recover from ()A.1|) the continuum 
WTI. Indeed, (|A.1|) can be renormalized to yield 

ZA{r){d;Alix)P\0)) = (A.2) 

= 2[Mo-M{r,Mo)]{P\x)P'm , + 0(a) 

(r,Mo) 



The two sides of eq. ()A.2|) possess a well defined continuum limit, as soon 
as they are multiplied by a suitable (logarithmically divergent) renormaliza- 
tion constant (i.e. the renormalization constant of the pseudo-scalar quark 
density), if the quantities Za and M are set to their appropriate values. 
Comparing the expression (|A.2|) with the form of the continuum WTI 

(a^i^(x)P^(O)) = 2m,{P\x)P''m , (A.3) 

one concludes that Za is the renormalization constant of the lattice axial 
current and the critical mass is the value of Mq which solves the equation 
Mq = M(r, Mq). We note incidentally that from the transformation prop- 
erties of the various terms in eqs. ()A.1|1 and (|A.2|1 under TZ^, one can infer 
the relations 

ZA{r) = ZA{-r) , (A.4) 
M(r, Mo) = -M{-r, -Mq) . (A.5) 

Eqs. ()A.4|) and ()A.5|) confirm the result that the renormalization constant 
of a multiplicative renormalizable operator is generally even under r — > — r 
(see sect. 13. 3|) . while the critical mass is odd (eq. (|2.7|) ). 

The relations between the renormalized continuum quantities indicated 
with the symbol ' in eq. (|A.3|) and the corresponding bare lattice ones are 

= ZaAI (A.6) 
= ZpP^ (A. 7) 

rhq = ZmZp^{Mo - Mcr) = Z'^Ulq , (A.8) 



It is, however, always understood that derivatives have to be thought as acting outside. 
As in the main text, we indicate nowhere the dependence upon the bare coupling constant, 
go, unless necessary to remove possible ambiguities. 
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where Z50 is the multipUcative renormahzation constant of the (subtracted) 
isosinglet scalar quark density and 



Zmir) = 1 
The last equality in eq. 



dM{r, Mo) 



Mo=M„(r) 



dMo 

follows from the identity 
_Zp_ 

Z5O 



(A.9) 



(A.IO) 



proved in ref. [H^ . 

As we argued in the text, in order to discuss the issue of improvement 
one has to make reference to the notion of effective theory introduced by 
Symanzik and use the related Symanzik expansion of lattice correlators in 
terms of correlators of the continuum theory. Recalling the analysis carried 
out in ref. [S], one gets up to 0(a) corrections the explicit expansions (x / 0) 



'■PP 



{r,mq) 



cont 



(mq) 



+ 



+avZ{r){P'{x)P''{0) J d'yCswiy)) 
+amqr^PP{r){P\x)P\<d) j (fyFF{y)) 
+amlri^P{r){P\x)P\0) / AV^^(y)) 



+ 



cont 

(mq) 

cont 



+ 



+ 0(a2 



(A.ll) 



{d,A'^{x)P'm 



'-AP 



{r,mq) 



cont 



+a4^ir){d,d,P'>ix)P''m 



(ruq) 



+ 



cont 



(mq) 



+ 



cont 



amgV^^ir){d,Alix)P\0) / d^yFF(y)) + 



"q 'I'tpip y I \^i^^pi. 

where FF is a shorthand for Tr[F^,yF^,y] and 



(rnq) 



+ 



(A.12) 



(A.13) 
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To avoid trivial 0(a) effects we have replaced the backward lattice derivative 
of eq. (|A.1() with the symmetric one, = {d^ + d*)/2. As in the text, by 
slightly changing the notation employed in eqs. HA.1|) and HA.2|) . we have 
indicated the parameters that specify the fermionic action using, besides r, 
the excess quark mass, ruq = Mq — M^, instead of Mq. In order to get 
the Symanzik expansions in the form given above, use has been made of the 
equations of motion and the continuum WTI (|A.3|) in the r.h.s. of eqs. (|A.11|) 
and (|Jl2]) . 

As explained in the text (sect.|31), the coefficients (, ^ and ry are all finite 
functions of r and and are necessary to match the lattice correlators 
on the l.h.s. to the continuum Green functions on the r.h.s. of eqs. HA.11|1 
and (|Jl2]) . There is an obvious correspondence between the appearance 
of the terms with the coefficients ^ and rj in front in the above Symanzik 
expansions and the need to introduce the coefficients cswj ^mj 

bq and ca, 

bA, bp to improve the action and the operators A^^, P^. 

Even in the limit of vanishing lattice spacing, eqs. HA.11|) and HA.12|) 
have a logarithmically divergent a-dependence which can be easily disposed 
of by multiplying the two sides of the equations by the appropriate renor- 
malization constants. 

We remark that in the expansions ()A.11|) and (|A.12|) . besides correlators 
of genuine multi-local operators ("fully localized operators" in the follow- 
ing), one has also to include terms where the multi-local operator under 
consideration (be it P^{x)P''{0) or 9^^|^(x)P''(0)) is inserted together with 
space-time integrated local densities. The former are related to the improve- 
ment of the various local factors making up the operator in the l.h.s., while 
the latter owe their origin from the need to improve the action [S]. Inte- 
grated terms require a little discussion, as the integration variable, y, will 
inevitably hit the points where the bi-local (in the present instance) opera- 
tor in consideration is localized (i.e. x and 0). To analyze the situation one 
may imagine to split the space-time integral in three different regions: one 
around x, a second one around plus the rest. Terms coming from inte- 
gration regions concentrated around the points x and can be shown, by 
an analysis based on the O.P.E., to lead to 0(a) contributions of the same 
form as those of the fully localized multi-local operators already present in 
the expansion. They only have the effect of redefining the corresponding 
Symanzik coefficients. The integration over the remaining portion of the 
space-time is finite and gives rise to genuine 0(a) contributions. 

Starting from the expansions HA.11|) and (|A.12|) , there are two equivalent 
ways to construct 0(a) improved lattice quantities, which we now want to 
make fully explicit: the Wilson average and the mass average method. 
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The Wilson average method 

Under the change of variables , we obtain the lattice formulae (valid for 
any value of a) 



{P\x)P\0)) 



{P\x)P\0)) 



{d,Ai{x)p'm 



i-r-Mo) 



(r,Mo) 



-{d,A'^{x)P\Q)) 



(A.14) 
(A.15) 



Using (|2.7j) . these relations can be rewritten, with the same little abuse of 
notation as before, in the more expressive form 



{P\x)P\Q)) 



{r,mq) 



{p\x)p'm 



{d,Aiix)P\0)) 



i-r-mq) 
bf 



{r,mq) 



{d,A'^{x)P\0)) 



{~r,—mq) 



(A.16) 
(A.17) 



To write down the equations that are obtained from the expansions HA.11|) 
and ()A.12|) when both r and mg change sign we refer to the analysis pre- 
sented after eq. H3.4() in sect.|31 

For completeness we list below the obvious parities of the continuum 
correlators appearing in the r.h.s. of eqs. HA.11|) and HA.12|) under nig 
—niq, which are implied by the transformation 



{p\x)p'm 

{p'{x)p'{o) J A/:sw(y)) 

{P'{x)P''{0) [ d'yFF{y)) 



{P''{x)P''{0) <fymv)) 



{d,A\xx)p'm 

{d,d,P\x)P\Q)) 
{d^Al{x)P\Q) J d'^yCswiy)) 
{d,Al{x)P\Q) I d'yFF{y)) 



cont 
(ruq) 


even 


(A.18) 


cont 
{rriq) 


odd 


(A.19) 


cont 
(m,) 


even 


(A.20) 


cont 


odd 


(A.21) 


cont 

{rUq) 


odd 


(A.22) 


cont 
(mq) 


even 


(A.23) 


cont 
{mq) 


even 


(A.24) 


cont 
{mq) 


odd 


(A.25) 


cont 
{mq) 


even 


(A.26) 
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Using equations from ()A.18|) to (|A.21|I . the identity l|A.16|) leads to the 
relations 



PPf 



(A.27) 



am,[epUr) + epU-r)]{P\x)P\Q)) 



-pp 



cont 



(mq) 



+ 



cont 



+ 



+amg[r^^kr)+7^^^i-r)]{P''{x)P''{0) I <fyFF{y)) 



^PP 



{rriq) 
cont 



+am'[rj?^{r) + r,Z{-r)]{P%x)P''{0) I d^y^^) 



^PP 



(iriq) 
cont 

{rriq) 



+ 



(A.28) 



from terms of 0(0*^) and 0(a), respectively. Similarly, using equations from 
(|A.22|) to (|A.26|) . the identity (|A.17|) leads to the relations 

Ci^r) = Cik-r) , (A.29) 



am, idkr) + Cii^{-r)]{dMx)P\0)) 



'-AP 



cont 



+a [rj^^r] + r^f,^i-r)]{d,d,P''{x)P\0)) 



^AP 



(niq) 
cont 



+ 



(mq) 



+ 



+a[4^{r)+r,i^{-r)]{dMx)P\Q) I d^y Cs^{y)) 



„AP 



cont 



+ 



+am, [vVF{r)+4F{-rW,Al{x)P\^) J d^yFF{y)) 



(rriq) 
cont 



+aml [rj^P(r)+rj^:;{-r)]{d,A''^{x)P'>{0) I d'y^ij{y)) 



-AP 



cont 



+ 



0. (A.30) 



From the above four equations one immediately comes to the conclusion 
that the average of correlators computed with lattice actions having Wilson 
terms of opposite signs and the same value of = Mq — Mcr are free of 
0(a) discretization effects. In formulae one gets, in fact 



1 r 



{P\x)P''{0)) 



{r,mq) 



+ {P\x)P\<d)) 



i-r,mq) 



C^kr){P\x)P\0)) +0(a2) 



(A.31) 



{d,A';,{x)P\0)) 



(r,mq) 



+ {d,Ai{x)p'm 



{-r,mq) 



Cikr){dM)P%0));-'+O{a'). 



(mq) 



(A.32) 
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Though unnecessary to obtain the improved formulae ()A.31|) and ()A.32|) , we 
observe that the vanishing of the sum of terms in eqs. HA.28|) and HA.30|) for 
arbitrary values of x implies that all the Symanzik coefficients rj and ^ are 
odd functions of r. 



The mass average method 

An equivalent way of getting 0(a) improved lattice data from standard 
Wilson simulations consists in comparing the two sets of formulae that are 
obtained by only changing sign to mq in eqs. (fOTTl and (f02ll . Taking 
into account the 75-chirality properties of the continuum correlators listed 
in eqs. ()A.18|) to ()A.21|) and (|A.22|) to (|A.26|) . it is immediate to verify that 
lattice correlators computed with fixed value of r and opposite excess quark 
mass, mq = Mq — Mcr, if linearly combined with a relative sign equal to 
the 7^5-parity of the multi-local lattice operator under consideration, yield 
quantities that are free of 0(a) discretization effects. In the case at hand 
one gets 



{P\x)P\Q)) 



{r,mq) 



C^^ir){P\x)P\0)) 



(mq) 



0(a2) 



(A.33) 



{d,A''{x)P\0)) 



(r,mq) 



{d,A'i{x)p'm 



{r-mq) 



afWMW^'(O)), ^ + 0(a^). 



(A.34) 



Eqs. (|A.33|) and (|A.34|) are perfectly equivalent to eqs. (|A.31|) and ()A.32|) . 

respectively, as it is immediately seen by using the lattice identities 



{P\x)P\{))) 



{r-mq) 



{P\x)P\0)) 



{-r,mq) 



{d,A'{x)p'm 



{r-mq) 



{d,Ai{x)p'm 



{-r,mq) 



(A.35) 
(A.36) 



which are direct consequences of ()A.16|) and (|A.17|) . 



Appendix B - Reflection positivity and transfer matrix 

In the continuum Euclidean theory the Wightman reconstruction theorem 
is based on the axiom of reflection positivity (HH, i.e. on the existence of 
an anti-linear operator, 0, which transforms an arbitrary functional of the 
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fields, F^, with support in a region, R, of the half-space t > into a func- 
tional, ©-F+, with support in the reflected region, 9R (with the definition 
9x = 0{x,t) = (x, —t)), such that 

{F+eF+)>0. (B.l) 

This condition implies the existence of a positive definite, self-adjoint Hamil- 
tonian. For a gauge theory the theorem is proved by choosing the gauge 
Aq = (temporal gauge). 

A similar theorem holds on the lattice [351 1361 07| . However, one has 
to distinguish between link and site reflection positivity. The associated 
reflection operations act differently on lattice points: the former, O^, refers 
to reflection across a hyperplane which cuts time links in half (e.g. the links 
(x, 0) — (x, a)), while the latter, Og, to reflection across a time hyperplane 
which passes through sites (e.g. the sites (x, 0)). Explicitly in the above 
examples one has 

9i{x,t) = {x,-t + a) e,(x,t) = (x,-t). (B.2) 

The validity of both site reflection positivity (eq. ()B.1|) with Q = @s) and 
link reflection positivity (eq. ()B.1|) with = 0^) is a sufficient condition for 
the existence of an Hermitean positive lattice (single-link) transfer matrix. 
Viceversa, if such a transfer matrix exists, one can prove that both types 
of reflection positivity hold. If only one kind of reflection positivity is valid 
(e.g. link reflection positivity, which typically holds under weaker conditions 
on the action parameters than site reflection positivity) one can still prove 
the existence of a positive transfer matrix acting over a time translation of 
two lattice sites (to be interpreted as the square of the single-link transfer 
matrix) . Naturally in the limit a ^ the condition of link (or site) reflection 
positivity alone is perfectly sufficient to extract continuum physics from 
lattice data. 

As in the continuum theory all these statements are proved in the tem- 
poral gauge, which in the lattice corresponds to set all the link variables 
associated to temporal links equal to the identity element of the gauge 
group (except at one set of links, which can be conveniently chosen to be 
(x,0)-(x,a)). 

For completeness we report the explicit action of the anti-linear reflection 
operation on the fundamental fields of the theory. Qg and Qi act on ipix), 
ip{x) and Ufj_{x) in the same way, the only difference is how the arguments 
they depend on get reflected. With the subscript s or £ understood one 
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has 



e[Uo{x)] = u^{ex-a6), e[Uk{x)] = u^{ex), k = 1,2,3 (b.s) 

e[^{x)] = V5(^x)7o , eii^ix)] = 7oV(^x) (B.4) 
and on monomials of fermionic fields 

G[/(C/)V(xi) . . . ^P{xn)] = ne[Umi}{xn)] . . . Qmxi)] , (B.S) 
f{U) being a functional of link variables only. As usual the definition 

^\ex) = e[$(x)] (B.6) 

is employed in this work. 

Site reflection positivity for Wilson fermions 

In ref. |HH| it was proved that site reflection positivity holds for Wilson 
fermions, provided 

\8r + 2aMo\>6, r = l. (B.7) 

The explicit construction of the transfer matrix of the lattice theory was 
carried out in ref. 1511 . 

The proof of site reflection positivity given in |3H] for r = 1 can be 
extended to r = — 1 by noting that commutes with TZ^. In fact, let 
Fj^{'il),ip,U) be a functional of the field variables defined at positive Eu- 
clidean times. Consider the expectation value {F^ QgF^) . The change of 
integration variables induced by the transformation leads to the rela- 
tion 



{F+{^P,i;,U) QsF+{ij,ij,U)) 



(-l,Mo) 



= (F+(-V;75,75V',C^)0si^+M75,75V',C^)) (B.S) 

(1,-Mo) 

Thanks to reflection positivity which holds for r = 1, one concludes that the 
r.h.s. of eq. (|B.8|) is non-negative, provided (eq. ()B.7|) with Mq — Mq) 

|8-2aMo|>6. (B.9) 

The inequality (|B.9() is precisely the condition HB.7|) , only written for r = — 1 . 
This shows that Wilson fermions obey site reflection positivity for both r = 1 
and r = — 1 under the condition |8r -|- 2aMo\ > 6,r = ±1. Link reflection 
positivity was discussed in sect. 12.31 

^■^In the two argument notation U^{x) = U{x,x + ajl) the action of <d reads Q[U{x, y)] = 
U*{9x, 6y), y = X + ap,. 
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Site reflection positivity for tm-LQCD 



The situation for tm-LQCD is similar to that of the standard Wilson theory 
summarized above. The validity of site reflection positivity was indirectly 
proved in [SJ, where along the lines of the work of ref. |54j the transfer 
matrix for the action ()4.15|) was constructed under the condition (uj = uJm + 

UJr) 

\8r + 2aMcr{r) + 2amg COS uj\ > 6, r = l. (B.IO) 

From this result the validity of site reflection positivity for the form ()4.18|) 
of the tm-LQCD action can be proved by noting that Qg commutes with 
the chiral fermion field transformation H4.3|) . One finds, in fact, under the 
condition ()B.10|) 

(LUr=Ul, UJm,=0) 



(F+ (^ph , Vph , t/) e, F+ ( V'ph , V'ph , f/) > 



(r,m,) 



with 



= (F+(V;,V,C/)esF+(^,V,^7)) J"" "^>0 (B.ll) 

(r,mg) 

= e-*"^^"^/Vph = V^phe-'"^^"^/^ (B.12) 

where the last inequality follows from the validity of site reflection positivity 
for the action (|4.15|) . 



Appendix C - Proof of eq. (13. lip 

In this Appendix we want to show that the 75-chiralities of the operators ap- 
pearing in the Symanzik expansion (|3.3j) , or in the analogous expansion (|5.1|1 
if tm-LQCD is employed, satisfy the equation 

Pn, [O] + Ptz, [Oe] +n, = l mod (2) . (C.l) 

Eq. (|C.1|) is an immediate consequence of the dimensional relation 1)3. 4() and 
the formula 

Pn, [O] + d[0] = Pn, [O,] + d[0,] mod (2) , (C.2) 
which we prove below. From these equations one gets, in fact, successively 

Pn, [O] + Pn, [Oi] +n, = Pn, [O] + d[0] + Pn, [O,] - d[0,] + 1 = 
= Pn^ [O] + d[0] - Pn, [Oe] - d[Oe] + 1 mod (2) = 1 mod (2) . (C.3) 

In order to prove eq. (|(y.2|) we note that with the definition 

A,[0] = PnM + d[0] (C.4) 
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eq. HC.2|1 becomes 

A5[0]-A5[0,] =0mod(2). (C.5) 
Eq. (|C.5|) is a direct consequence of the following symmetry considerations. 

1) - The Wilson and tm-LQCD actions (eqs. 1)2. 2|) and (|4.18|) . respectively), 
and a fortiori the formal continuum action, are invariant under the combined 
transformation TZ^ x V^, where '^^ 

U^{x) Uli-x-afi) 
Vd : < ip{x) e3*V2^(-x) (C.6) 

The invariance of the pure gauge part of the action is obvious. For the 
fermionic part invariance follows by first noting that under T>d (recall the 
definitions ((231)) 

^P{x)J2J,{V, + V;)^P{x) ^ ^p{-x)J2l,{V^ + V;)i;i-x) , 

Hx)Y.v;v^^p{x) ^ -^{-x)J2v;v^^p{-x) , 

■ip{x)ip{x) — > — 'ip{—x)ip{—x) . (C7) 

The minus sign in front of the chirally non-invariant terms is compensated by 
the transformation TZ^. Finally the space-time arguments are brought back 
to their initial values by the trivial change of space-time integration vari- 
ables —x — > X in the action. In the formal continuum theory an analogous 
argument goes through. In fact, the continuum version of the transforma- 
tion Vd, besides changing the sign of all the space-time coordinates, has 
simply the effect of multiplying each local term Ci in the action density by 
the factor (— l)'^'^'!, where d[Ci] is the naive dimension of £j. 

2) - Performing the change of functional integration variables induced by 
the transformation TZ^ x P^, one can prove that the lattice correlators of the 
type we are interested in obey the formula 



(^) 



{0{xi,X2,...,Xn))]''' = (-1)^^[^](0(-Xi,-X2,...,-X„)) 7' ^ (C.8) 



{r,mq) 



■^^In this Appendix, for notational uniformity, we will always denote the quark fields by 
■0 and ^ (as in eq. I2.2ll 'l. even when we refer to the tm-LQCD action (14.1811 . 
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with w = in the case of the standard Wilson action 1)2.21) . The rela- 
tion (|C8|) directly follows from our definition of A5 (eq. (|C4() ) and the way 
Vd (eq. (iniSl)) acts on the fundamental fields in O. 

Similarly, all continuum operators contributing to the Symanzik expan- 
sions (|3.3|) or (|5.1|) of the expectation value of 0{xi,X2, • • • , Xn) must obey 
eq. HC.8|) with a value of A5 equal (mod (2)) to that of O. In other words all 
terms in the Symanzik expansion of O must have the same (—1)^^ -parity 
as O. This is just our thesis, eq. (|C.5|) . 



Appendix D - Discretization effects in the determination of 

M„ 

In this Appendix we want to show that 0(o) discretization effects in the 
determination of Mcr(r) do not jeopardize the proof of 0(a) improvement 
of averaged correlators given in the text. We give the proof in the case of 
standard Wilson fermions. The generalization to tm-LQCD is straightfor- 
ward. 



The problem 

In PT the analytic dependence of lattice correlators on a is completely and 
explicitly computable. Thus an unambiguous expression of Mcr(r) can be 
given to any order in gQ. The situation is different non-perturbatively. An 
estimate of Mcrir) can be obtained by using, for instance, PCAC or the 
observation that (in the presence of S^SB) the pion mass is expected to 
vanishes at the chiral point (modulo cutoff effects, see sect. 12. 2|) . If the the 
pion mass is employed to determine Mcr(r'), one has to proceed by extracting 
from the large t behaviour of the two-point correlator -^"^ 

G^(t;r,Mo) = a3^(P^(x,t)P^(0)) (D.l) 

X ir,Mo) 

a lattice measurement of the pion mass, m.,^{r, Mq). A sensible estimate of 
the critical mass, which will be referred to as M^^ (r) and obeys 

M^Hr) = -M^H-r), (D.2) 

can thus be obtained as discussed in sect. 12.21 

■^^In this appendix we label lattice expectation values by using, besides r, the bare quark 
mass. As in the rest of the paper, the ubiquitous Qq dependence is understood. 
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If PCAC (see e.g. eq. (|A.2|) ) is employed to measure Mcr(r), one would 
get another, equally good, non-perturbative estimate of it, say M^^, given 
by the solution of the equation 

M(r, (r)) - M^^ (r) = . (D.3) 

The new estimate is again odd in r (see eq. ()A.5|) ). 

For the sake of this argument we may summarize the situation by for- 
mally writing for any such estimate 

M^,{r) = M,,{r)+6'{r), (D.4) 

with 

S%r) = -5%-r) = 0{a) . (D.5) 

The difference Mcr{r) — M^^r) is obviously odd under r — > — r, because 
Mcr{r) as well as anyone of its estimates is odd in r (see sect. EJ. 

The problem we want to address is whether 0(a) ambiguities in the 
knowledge of the critical mass can invalidate the argument for 0(a) im- 
provement developed in sect. 01 leading to eqs. (|3.9|) and (|3.12|) (or to the 
similar ones valid in the case of tm-LQCD). We show below that eqs. (|D.4|) 
and ()D.5|) are indeed all is needed for eqs. (|.S.9|) and (|3.12|) to be valid. 

This result should not come as a surprise, since effectively we are here 
merely claiming that all the 0(a) cutoff effects, which generically are present 
before averaging, get canceled in VK^'s of correlators and derived quantities. 
Among these cutoff effects there are, indeed, also those arising from 0(a) 
errors on Mcr(r). They are not in any sense special and actually stem from 
the presence of 0(a) cutoff effects in the correlators employed to estimate 
the value of the critical mass. 



The proof 

Once the values of r (> 0), Qq and the bare lattice mass Mq = Mq are 
chosen, the Symanzik expansion of the expectation value of the operator O 
takes the form (see eq. ()3.3() in sect. ^ 



{O) 



(r,Af+) 



q 

cont 
«) 



+ 0(a^ 



«) 



+ 



where 



(D.6) 
(D.7) 



60 



is the continuum mass at the subtraction point a~^. is defined in terms 
of the (unknown) value of the critical mass, Mcr(r). 

We want to average this expectation value with the similar expression 
obtained by changing sign to r and adjusting Mq to a suitable value, Mq — > 
, such that all 0(a) terms cancel in the sum. We have shown in sect. El 
that this is indeed possible with the choice 

Mq- - Mcr(-r) = M+ - M^,{r) = m+ . (D.8) 

The "practical" problem with this prescription is that, as we just recalled, 
one can only get an "estimate" , M^j. (r) , of the critical mass with the prop- 
erties summarized in eqs. ()D.4p and (|D.5|) . 

Now the question is: given Mq , which value of the bare lattice mass do 
we have to take in practice, when we go over to — r, in order to get rid of 
0(a) terms? 

Clearly the only possible answer (and, as we shall see, the right one) is 
to consider the expectation value of O computed with —r and a bare mass, 
M'o , given by 

M'o" ^ M+ - M'Jr) + MtA-r) = M+ - 2M^,{r) , (D.9) 

where the last equality follows from the r-parity properties of M^^r). The 
Symanzik expansion of {0)\(^_^. takes the form 



(O), ^,,.=[CS{-r) + am-£g{-r)]{0) 



cont 



_ — , ^ cont „ 

+a J2i^-rvSA-r){0e) ^ _^ + 0{a') , (D.IO) 
with the definition 

m- = MV - Mer(-r) . (D.ll) 

For future use we note the chain of equalities 

m- = M'o - Af„(-r) = M+ - 2M^,{r) - M,,{-r) = 
= M+ - 2M'Jr) + Mcr(r) = M+ - M,,{r) - 25^ {r) = 
= m+ - 26^{r) , (D.12) 



where eqs. ()D.9|1 . ()D.4|) and (|D.5|) have been employed. 

Clearly the 0(a) mismatch between Mq (eq. ()D.8|) ) and M'g (eq. ()D.9|) 1. 
as well as that between and m~ matters only for the first term of the 
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expansions HD.6|) and HD.10|) . because all the other terms have an explicit 
factor of a in front of them. Thus summing HD.6() and (|D.10|) all 0(a) 
contributions that are multiplied by the (r odd) Symanzik coefficients ^ and 
rj cancel and we get (see eq. p.Tf) ) 



(r,M+) 



+ {o) 



C8ir) (O) + (O) 



cont 



«) 



cont 



(m, ) 



+ 0(0^). (D.13) 



At this point, recalling the expressions ()D.7|) of and HD.12|) of , we 
can put the r.h.s. of the previous equation in the form 



C8ir) (O) + (O) 



cont 



«) 



cont 



(m, ) 



C8ir) (0) 



cont 



(M+-Mc,(r)) 



+ {o) 



+ O(a^) = 

cont 

(M+-Afcr(r-)-25<=(r)) 



+ 0(a2). (D.14) 



In this expression the two continuum correlators are evaluated at "slightly" 
different values of the excess quark mass, namely Mq — Mcr(r) and Mq — 
Mcr{r) — 25^ (r). Neglecting O^a^) terms, they can be immediately brought 
to the common value of the mass parameter 



rUg = M+ - M,,{r) - 5%r) = M+ - M^,{r) 



(D.15) 



by a simple Taylor expansion of the continuum correlators in powers of 5^{r). 
The terms with the mass derivative of continuum correlators cancel out and 
one gets 



1 



(O) 



{r,mq+Mi^(r)) 



+ {0) 



{-r,m,-Mi,{r)) 



Coir){0) ; ; + 0(a^), (D.16) 



cont 



(mq) 



where only well defined, numerically accessible inputs, like M^^[r) and the 
excess quark mass, niq, appear. 

We end by noticing that, had we used the excess quark mass to label 
lattice correlators (as we did in the text after sect. the formula (|D.16|) 
would have actually read 



1 



(r,mq) 



+ (0> 



(—r,mq) 



cont 



C^(r)(0) +0(a^) 

[rriq) 



(D.17) 



The result of this elaborated discussion is that 0(a) improved correlators 
are obtained using the formula ()D.17|) - which is identical to the W A pre- 
scription ()3.9|) - with niq given by eq. ()D.15|) . 
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Appendix E - Improvement of the baryonic spectrum 

In this Appendix we want to give a quick argument to show that the pro- 
cedure for 0(a) improvement proposed in this paper (eq. (|3.26|) and its 
tm-LQCD counterpart) works also for the baryonic sector of the hadronic 
spectrum. For concreteness we wih concentrate on the spin 1/2 and 3/2 
states that can be formed with an SU{2){ flavour doublet of quarks. 

First of all we note that all baryonic energies can be extracted from 
spin averaged two-point correlators of the type 1)3. 18(1 and ()5.1U|) . for the 
untwisted and twisted Wilson formulations, respectively. The (bare) inter- 
polating operators for spin 1/2 and 3/2 baryonic states are of the form 

iHji^) = [^^Jli^LM]i^), (E.l) 

where / and a are free flavour and Dirac indices, I represent (contracted) 
color indices and is a bosonic operator with definite transformation 
properties under TZ^. For instance, for proton {f = 1, = u) or neutron 
(/ = 2, ip'^ = d) states we can take 

M'ln = i^jl5r24^i (E.2) 

with 'ip'^ denoting the charge conjugated of "ip. Similarly, the (bare) interpo- 
lating operator for the spin 3/2 baryon A"*"^ is of the form (|E.1|) with / = 1 
and 

-^Ifc = ^(1 + ^3)7fc i^i, k = 1,2,3. (E.3) 

In all cases the bi-local operator ^h{x)^\{0) (with all Dirac indices con- 
tracted) has definite (positive) T^s-parity. Since this property is preserved 
also in presence of mixing, we can apply the considerations of sect. 13.21 (or 
the analogous ones valid for tm-LQCD) and conclude that also baryonic 
masses are improved by the usual WA procedure. 

Moreover, in the tm-LQCD formulation with action 1)4. 18() the masses of 
the baryonic states HE.1|) are even functions of the twisting angle lo, hence 
automatically free of 0(a) corrections if a; = ib7r/2. To prove this result it is 
sufficient to check that the argument presented in sect. 15. 3 ITI for the w-parity 
of masses applies to the states under consideration here. This is actually so, 
because the interpolating operators (|E.1|) are such that the spin averaged 
product $/i(x)<I>|^(0) has obviously a definite positive formal parity (eq. ()5.5|) ) 
and this property is preserved under renormalization. Incidentally, itself 
could be chosen to have definite formal parity, if one decides to multiply the 
operator in the r.h.s. of ()E.1|) by the projectors (1 it 7o)/2. 
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Appendix F - Proof of eq. (j5.18j) 



A detailed proof of eq. (|5.18|) requires repeating in the case of tm-LQCD the 
hne of arguments that in a parity invariant theory leads from the way V acts 
on the elementary fields entering the Lagrangian density (see eq. ()5.4|) ') to the 
notion of parity for the eigenstates of the Hamiltonian. Observing that in the 
case at hand (action 1)4. IBI) ) the spurionic symmetry V^^ = V x {uj ^ —to) 
((7?sp)2 _ -j^^ replaces V, one can proceed as follows. 

1) To composite bare operators formal parity labels are attributed ac- 
cording to the naive prescription induced by the field transformation (|5.4j) . 

2) A basis of m.r. operators can always be constructed, consisting of 
operators endowed with a definite parity label. In fact, to comply with the 
symmetry V^^ the mixing pattern will have to be such that (bare) operators 
with opposite formal parity get mixed with coefficients whose relative uj- 
parity is negative. 

This procedure is equivalent to assign a parity label via eq. (|5.5|) . If 
one wishes to employ this method, one should start by assigning parity to 
local operators (e.g. quark bilinears) or simple products of local operators 
with non-zero (for suitable choices of the space-time arguments) vacuum 
expectation value and then move on to progressively more complicated cases. 

The situation here is exactly analogous to what happens with the spu- 
rionic symmetry TZ^^ , where eq. 1)3. 2|) is exploited to attribute a definite 
75-chirality, P-jir^lO], actually the naive one, to any (gauge invariant) m.r. 
composite operator, O. 

3) For the vacuum state we have in general 



v\n) 



{r,mq) 



■r]^{uj)\^l) 



{r,mq) 



(F.l) 



with rjQ{tLi)riQ{—uj) = 1 (from "P^ = l |S3]) and |r7Q(u;)p = 1 (from the 
vacuum normalization condition). In view of these results one can always 
decide to simply take the phase of the vacuum state such that i^niuj) = 1 for 
any value of uj. If one does not want to do so, one will have to keep track of 
the vacuum phase in all the forthcoming formulae. 

4) A basis for the states of the Fock space is constructed by applying 
creation operators to the vacuum (cyclic vector). The elements of the basis 
are highly non-trivial and schematically can be written in the form 



k; a 



(■r,mq) 



{r,mq) 



(F.2) 
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where, according to 1) and 2), the creation operators a|j^(k) (a is an index 
which labels the states with the same unbroken quantum numbers h and 
three-momenta k) will satisfy the relation 



^«L(k)^ = ^/^"«L(-k). 



(F.3) 



where r]ha (= ±1) is by construction an to- and k-independent label. 

5) The eigenstates of the transfer matrix T{uj,r,mq) (see eq. (|5.13() 'l will 
be expressed as linear combinations of the states ()F.2|) of the kind 



\h, n, k) 



I {r,mq) 



Cna {uj,r,mq)\ h, k; a) 1^^^ 



(r,mg) 



(r,mq) 



(F.4) 



with appropriate coefficients, Cna{'^^f,rnq), to ensure normalizability of the 
state and consistency with the spurionic invariance P^p. 

In particular, consistency with the spurionic symmetry V^^ implies the 
operator relation 



V 



\h, n, k)(/i, n, k| 



(r,inq) 



V 



n, — k) (/i, n, — k| 



(r,mg) 



(F.5) 



which is proved below in subsect. ll^'.ll This relation can only hold if the co- 
efficients CnQ{uj,r,mq) have (up to an irrelevant a-independent phase, which 
we set equal to unity a definite parity, Sna = il, under lo — > -co, i.e. 
only if they satisfy the relation 



,{uJ,r,mq) = SnaCna{-i^,r,mq) 



with 



Vh,r. 



an Q-independent label. 

6) For the states ()F.4|) the relation 

P|/i,n,k)j^^J^_^^ = r]h,n\h,n,-k] 



{r,mq) 



(F.6) 
(F.7) 

(F.8) 



is immediately seen to be a consequence of eqs. (|F.3() . (jF.Gf) and HF.7() . By 
construction (see eq. ()F.7|) ) rj^^n can only take the values ±1. The rela- 
tion (|F.8() is nothing but eq. ((^7TH|) of sect. 15.3.21 

^^This phase can always be readsorbed in a physically unconsequential redefinition of 



the overall phase of the state n,k)| 
section. 



(r,m,) • 



See also the comment at the end of this 
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Whether the parity label r]h^n is 1 or —1 can be decided numerically by 
the method discussed in sect. 15. 3. '51 The resulting value of % „ - at least for 
values of the bare parameters for which quantities of O(a^) and 0{a^) can be 
unambiguously distinguished - is to be identified with the parity quantum 
number of the state corresponding to |/i,n, k)(^ „j ) in the continuum limit. 
By construction rjh^n is independent of the three-momentum of the state, but 
does depend on the excitation level n because V is not an exact symmetry 
of tm-LQCD. There may exist, in fact, states with opposite parities but the 
same unbroken tm-LQCD quantum numbers. For instance in the vacuum 
sector one encounters the vacuum itself (with ijq = 1), the neutral pion 
(with r]^o = — 1), etc. 

We observe in conclusion that one could decide to separately alter the 
phases of the states \h,n,'k)^'^^^ and \h,n, —k)^,'^ with respect to what 
one gets from the above construction. If one does so, then track of this 
modification should be kept in all the relevant equations, in particular in 
eq. ()F.8|) . The only consequence of this choice is that all lattice formulae 
relating quantities with opposite values of cj, or with values of ui differing 
by vr, become unnaturally and unnecessarily cumbersome. 



F.l Proving eq. ()F.5|) 

Let us finally come to the proof of eq. (|F.5|) . For this purpose it is convenient 
to consider the general two-point correlator 

G$,s,(k,i;a;,r,m,) = a35]e*'^-(f,(x,t)4(0,0)) J"^ (F.9) 

where $/i and S/^ are two arbitrary local operators with formal parities 7/$^ 
and r/s^ and equal unnbroken quantum numbers h. Performing the change 
of variables induced by (|5.4|) , one obtains 

G^,^s,S'^,t;uj,r,mq) = T/$,^T/H,,G$^Sft(-k,t;-w,r,mg). (F.IO) 

By comparing the spectral representations of the two sides of eq. HF.10|) and 
using eq. (|5.15j) . we can get the set of relations 



{n\<^hiO,0)\h,n,k){h,n,k\El{0,0)\n) 



t .M^v H 



r?$,r/H,(^^|$h(0,0)|/i,n,-k)(/i,n,-k|4(0,0)|J7) (F.ll) 

(r,mqj 
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thanks to the leverage we have from the dependence of G^^h^ on t. Inserting 
factors "P^ = 1 on the l.h.s. of eq. (|1\11|1 yields 



If we now combine eq. (|i:'\l() with eq. ()1'.12|) and exploit the arbitrariness of 
the operators and H/^, we arrive at the relation ()F.5|) . 
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